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I. 


INTRODUCTION 


The development of the LANDSAT earth resources satellite has resulted 
in the availability of new information which has immediate application 
and utility to environmental planners and managers. However, in order 
for this new information to be absorbed into common pratice, efforts 
must be made to transfer it from the research sector to the user sector. 

Two important steps are required for this transfer process to 
occur. 

1) The applicability and accuracy of individual remote sensing 
techniques must be confirmed in the user arena. 

2) The economic practicality of remote sensing in augmenting 
or supplanting conventional data sources must be clearly 
demonstrated . 

Ecosystems, under contract to the NASA Marshall Space Flight Center, 
has focused upon such technology transfer in the area of hydrology and 
water quality. The ECHOS hydrologic model and the hydrologically oriented 
remote sensing analysis procedures developed for NAS8- 30539 have formed 
the basis of the transfer task. 

Specifically, the work included four tasks: 

Task 1 . Demonstrate the hydrologic model developed under the 
previous phase of this work and its corresponding use 
of LANDSAT imagery to selected users. Explore with 
the users the active transfer of remote sensing. 

Task 2 . Perform demonstrations of the value of the hydrologic 
model and/or of extracting ground cover information 
from LANDSAT data as may be requested by the more 
important users. 
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Task 3, Publish the material, or present it at selected 

symposia, and similar meetings as indicated by the users. 

Task 4 . Document user reactions, specific comments, and 
positive actions taken. Assess and document the 
nature and magnitude of potential impediments to 
technological transfer. Generate a case history 
of user transfer experience and guidelines as an 
aid in similar efforts by NASA. 

True technology transfer is roost beneficially accomplished when the 
transfer encompasses the overall problem as contrasted to problem support 
areas. The users are mainly interested in the overall results that accrue 
from the use of remote sensing and secondarily in the specifics of remote 
sensing, per se. However, once the utility of the remote sensing 
application is proven, the diffusion of the remote sensing techniques 
within the user's organization rapidly occurs. Consequently, this 
project sought to address end-to-end solutions of important problems 
where remote sensing would play a key role. 

The major project chosen for the technology transfer demonstration 
was the estimation of non-point source water pollutant loading from 
storm runoff. The need for such data is national in scope and of 
immediate importance as a consequence of the federal "208" water quality 
program, which is described in Section III. LANDSAT derived land 
cover assessments and the hydrologic parameters therefrom is the key 
input to non-point source load analysis. 

Althoug?! "208" is a national program, a local focus exists in the 
Baltimore Regional Planning Council (RPC). This body is charged with 
preparation of the water quality plan for the 6,123 square kilometer, 
six county, central Maryland region surronding Baltimore. Demonstration 



of the transfer of LANDSAT analysis techniques and the LANDSAT data 
oriented ECHOS model at the RPC has both local value and national 
applicability . 

The objective of the project was to determine hydrologic response 
and consequent pollutant loads from important watersheds using remotely 
sensed land cover. In joint consultation with the RPC, ECOsystems was 
asked to apply LANDSAT techniques to the estimation of yearly pollutant 
loads for one of its three designated prototypical test areas, the 
Magothy Basin in Anne Arundel County, Maryland. 

The progress made during the past year has utilized the research 
and development previously performed for NAS8- 30539. The pertinent 
results of that contract are presented in the following section. 
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II. SUMMARY OF PREVIOUS WORK 

From February, 19 74 through February, 19 77, ECOsys terns has been 
performing research directed toward the development and validation of a 
LANDSAT data oriented hydrologic computer model under MSFC contract 
NAS 8- 350 39. The effort included the following investigations; 

1) Identification of those hydrologic phenomena (infiltration^ 
antecedent soil moisture, etc.) which ’’drive" peak runoff 
events, i.e., those parameters to which the watershed's 
outflow is most sensitive. 

2) The incorporation of these factors into a hydrologic 
computer model compatible to the maximum degree with inputs 
from remote sensing. 

3) The verification of the model for actual watersheds. 

4) The documentation of procedures for extracting hydrologic 
information from LANDSAT imagery. 

The ability to estimate streamflow from rainfall events is one of 
the primary requirements of water resources planners and engineers. It 
impacts the design and construction of waterworks and the plans for 
future uses for proximate land areas. 

The most useful tool available to water resources managers is the 
computer model designed to simulate the response of a watershed to a 
given rainfall input. For acceptable performance, these models require 
specific physiographic and hydrologic data about the watershed. Typical 
among model parameters are flow lengths, slope, surface friction, soil 
moisture, and channel characteristics. Current data is either unavailable 
or expensive to obtain for many basins; therefore, it is advantageous 
to exploit the low cost and multi— temporal features of satellite remote 
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sensing to provide more con^lete information as well as periodic up- 
dates, Many of the parameters required are surface or surface-infer- 
red phenomena and, therefore, are potentially- measurable from earth 
resources satellites. 

The result of contract NASS- 30539 was a hydrologic computer model 
oriented to the use of LANDSAT data. The model has been tested on well- 
instrumented watersheds with good results. A set of procedures was also 
derived for the extraction of necessary input data from LANDSAT images. 
These procedures involved the analysis of LANDSAT multispectral and 
multi- temporal images by means of an optical image analysis. The tech- 
niques were applied to land cover classification of typical Maryland 
watersheds. Average inventory errors were typically five percent of 
less for the important classes. These accuracies were completely satis- 
factory for derivation of hydrologic parameters such as Manning’s "n" 
and sub-surface abstractions. The computer model and the analytic 
techniques employed were described in detail in the final report 
(CR150235; 150236), 

While the utility of the techniques and procedures of rainfall to 
rmoff estimates were of definite interest to the important research 
users, it was obvious that a more practical technology transfer of re- 
mote sensing requires the extension of the techniques to the user com- 
munity. Examples of these user applications are: simulation of flow 

peak in a flood prediction application, or the use of the rainfall/ 
runoff to estimate pollutant loads. An important result from this 
analysis was the realization that LANDSAT data could have an important 
role in the nationwide EPA sponsored "208" program, discussed in the 

,:p\ following section. 

1 
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III. TECHI^OLOGY transfer to the BALTIMORE REGIONAL PLANNING COUNCIL 

The U.S. Environmental Protection Agency, under Public Law 92-500, 
Section 208, has adopted regulations mandating regional jurisdictions 
to develop areawide waste water management plans. The purpose of these 
programs is to ensure that the areas meet the projected 1983 and future 
water quality regulations. The plans will be used to guide water quality 
programs through the year 2000. National federal expenditures for "208** 
programs are presently set at $380 million over a two year time frame. 
Legislation is currently pending to project future funding requirements. 

Until recently, the major thrust of the EPA has been on point 
sources, i.e., those pollution sources which are confined to spatially 
limited areas. The "208" program has expanded the scope to include the 
non-point sources. In contrast to point sources, non-point sources result 
primarily from the impact of the hydrologic cycle with the land cover. 

Land cover and associated environmental conditions determine the type, 
form, concentration and temporal pattern of pollutants for each watershed. 

Quantification of non-point pollutant load requires: 1) current 

and frequent assessment of land cover over large geographic areas; 

2) the determination of the relationships between diverse land cover 
classes and the pollutants they produce; and 3) the determination 
of the impact of the hydrologic cycle upon the surface and groundwater 
transport of pollutants. The rapid classification of stable and transient 
land cover, therefore, is a necessary and valuable adjunct to the "208" 


program. 
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Approximately 9-10% of total federal "208" expenditures can be 
directly or indirectly related to land cover measurements and land 
cover/non-point source pollutant load correlation. Regional environ- 
mental planners with "208" responsibilities, therefore, are an impor- 
tant potential user class interested in LANDSAT derived hydrologic 
land cover and pollutant load data. In many cases, only the LANDSAT 
type data has the sufficiently good resolution and freshness necessary 
to accommodate the required loading estimates , 

Because of this potential, it was determined to undertake tech- 
nology transfer with a Section "208" water resources planning body. In 
the metropolitan Baltimore area, the "208" program is managed by the 
Baltimore Regional Planning Council (PJPC). The RPC is a multi-juris- 
dictional organization serving a section of the state of Maryland com- 
prised of the City of Baltimore, and the Maryland counties of Anne 

Arundel, Baltimore, Carroll, Harford, and Howard. The combined area of 

2 2 

the region is approximately 6123 km and is composed of 5841 km of land 

2 

surface and 282 km of surface water. The total area is equivalent to 
about one-fifth of a LANDSAT frame. Figure 1 shows the member juris- 
dictions and their proximity to the major population centers of Baltimore 
and Washington. 

The RPC is currently completing a two year, $2.5 million contract 
from the EPA. The object of this program is to develop pollution abate- 
ment plans. The effort is designed to provide the basis for a con- 
tinuing planning process within this multi-jurisdictional framework. 
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The RPC has an in-house staff of some 16 professionals sup- 
ported by 39 specially hired individuals in the separate city and 
county jurisdictions. The effort is advised on two levels: by pro- 

ject officers (6) from the jurisdictional public works or health of- 
fices and by the political heads (6) of the counties, assisted by the 
presidents of the county councils (6) . Figure 2 indicates the staff 
and advisory personnel and their affiliation to the RPC-208 project. 
The total project complement is 118, 60 of which are directly working 
on the project; the remaining 58 are public, local and federal govern- 
mental advisors. This is an important user audience for LANDSAT 
generated data products. 

Areawide Waste Treatment Management Planning under Public Law 
92-500 for the Baltimore area officially got underway in April, 1976. 
The Regional Planning Council, as the designated agency, was charged 
with the responsibility of preparing this areawide plan and for ad- 
ministering the grant funds. 

The Regional Planning Council contracted with each of the local 
jurisdictions to perform various elements of the project control pro- 
gram. The Regional Planning Council also contracted with the State of 
Maryland Water Resources Administration, to conduct non-point sampling 
program and to provide liaison between the state "303” Basin Plan work 
and the Regional Planning Council’s "208" program. In addition, a 
contract was awarded to the U.S, Soil Conservation Service to provide 
accurate estimates of soil erosion in the region. 
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The entire work program is being carried out under the direction 
of the Regional Planning Council's ”208'* staff. Each of the local juris- 
dictions participating in the areawide plan also has a ”208” staff 
headed by a project officer. There are a number of subcommittees com- 
prised of members from the RPC ”208” staff and the local jurisdiction’s 
208 staff. Each of these subcommittees meets on a regular basis at 
least once a month to develop the methodology for carrying out the 
work and monitoring the progress of the various elements of the re- 
spective section. 

A Public Advisory Committee (PAC) composed of interested citizens 
and a Technical Advisory Committee (TAC) composed of technical people 
provide input to the RPC and local jurisdiction’s "208” staff. A 
coordinating committee composed of elected officials from each jurisdic- 
tion sets the general policy and all program outputs are subject to 
the approval of this committee. 

The basic requirements of this initial phase in the ”208” planning 
process were: 1) to assess the present state of pollutant loads and 
pollution in the region identifying areas of criticality; 2) to de- 
velop techniques for quantifying the effect of various land uses on pol- 
lutant inputs to the receiving waters; 3) to predict the possible fu- 
ture growth of the region and its effect upon the key population and 
land use parameters which impact pollutant levels; and 4) to develop al- 
ternative structural and non-structural procedures to mitigate and 
control waterborne pollution to cost-effectively meet the EPA standards 
and guidelines. The development of an accurate , spatially-distributed 
data base is a key element to the successful operation of the ’'208" 
planning process . 



- 11 - 


The paucity of accurate data on land cover, water quality, and water 
quantity has made the technical aspects of the program extremely 
difficult. Conventional reliance on "black-box" models relating the 
water quantity and pollutant load outputs to precipitation inputs are 

neither reliable nor cost-effective. Simple "hand— type" computations made 
without taking into account the actual and current land cover can pro- 
duce an unacceptably high error level. However, due to the large 
areas under consideration and the high cost per drainage area for 
even the simpler models, greater reliance on an accurate measurement 
of surface cover and estimation of hydrologic properties therefrom is 
required. 

Under direction of the RPC and with MSFC concurrence, ECOsystems 
applied LANDSAT data for 1975 and the previously developed hydrologic 
model to the estimates of pollutant loads in the Magothy Basin, The 
Magothy is one of the three prototype basins selected by the RPC to develop 
techniques for region-wide pollutant load assessment. The following 
report details the application of the LANDSAT techniques to pollutant 
estimation in the Magothy. 
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IV. THE USE OF LANDSAT DATA IN THE ”208” PLANNING PROCESS 

Non-point source pollutant loads are dependent upon the hydrologic 
response of an area and the pollutants incorporated in and upon the surface 
of the area. Both the hydrologic response and pollutant accumulation 
are related to the nature and distribution of the land cover. 

Comprehensive data bases do not exist at present for specifying the 
exact cover of the land surface at a resolution sufficiently small to 
accommodate existing data on soils and sediment. Conventional land use 
maps do not have the resolution required to assess pollutant loads from 
individual small watersheds. The LANDSAT 70m pixel (1,1 acres) is 
sufficiently small, in most cases, to permit land cover accuracies 
commensurate with the estimates of pollutant loads. 

A basic land use inventory is maintained by the Regional Planning 
Council updated to 1973. Local zoning boundaries are also available, 
giving a composite representation of land use and zoning, albeit somewhat 
dated. This land use inventory does not comprehensively yield land 
surface cover. 

The land use data is currently being updated, reviewed, and 
validated by ”208" staffs and made available for use in ”208" planning. 

This inventory includes areas classified in the following categories: 
residential, commercial, industrial, institutional, open-space, 
agricultural, semi-public and military. 

In addition to the urban land use inventory, other data sources 
being utilized in the "208" planning process include current building 
permits, classified as residential (single family and multi-family) and 
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non-residential, suiranarized by small planning districts; an inventory 
of all land subdivisions in the region since 1970; and electric meter 
connection data, by 1/4 square mile grids, provided semi-annually. 

Regional "208" planning requires that the RPC and participating local 
governments upgrade their land use information data handling capability. 
Data files are now being created in digitized formats that enable computer 
mapping and graphic displays which expand land use analytic capability. 

Data files are being created for inventories of vacant land within 
planned sewer service areas, zoning restrictions, existing urban land uses, 
and other data on the region's land. 

Initial results of the application of computer classified LANDSAT 
data at RPC, while encouraging, did not yield accuracies commensurate 
with the program requirements for the small hydrologic basin levels. 

A land cover classification by computer was performed at RPC of 
the Rhode River Estuary, in Anne Arundel County, Maryland. Ground truth 
land cover data were supplied by the Smithsonian Institute Chesapeake 
Bay Center for Environmental Studies (CBCES) . The Chesapeake Bay 
Center's review provided land cover data for nine small sub-watersheds 
of the Rhode, ranging from 15 to 660 acres in size. The CBCES staff 
had complied their land cover data in five classes (forest, old fields, 
cropland, pasture, and residential land) . The fifty LANDSAT classes 
developed were aggregrated into combined classes for comparison. Dif- 
ferences between the CBCES and RPC results varied, depending on the 
class type, for each watershed. The weighted average difference be- 
tween the two data sources was 32 percent, not considered sufficiently 
accurate for the "208” application. 
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It was noted that the LANDSAT classifications used in the comparison 
were early April vintage, compiled prior to an e^diaustive reclassifica- 
tion of agricultural lands done" in cooperation with the Baltimore County 
Agricultural Extension Service and Soil District Conservationists. The 
CBCES data was from the late summer and fall of 1975 and represented 
photointerpretation amended by field checks. 

Following additional modification of the LANDSAT signatures, a few 
locational errors were corrected. The weighted average difference 
between the two data sources was revised to 20 percent using the revised 
signatures. Ihe twenty percent difference, however, was still seen as 
unacceptable in view of the modeling efforts planned and designed around 
the LANDSAT classifications. The required error limits for individual 
classes are of order 5-10% for sub-watershed areas of 1,000 hectares. 
These error bounds are required to maintain final estimates of pollutant 
loads in the range of 15 to 20% over the entire region. Subsequent 
adjustment of signature mixes brought the computer classifications within 
acceptable bounds of approximately 8%. 

Figure 3 indicates classification error bounds and the performance 
of the RPC computer-only classification schemes. In the individual 
sub-basins the classes considered were: land in crops, pastures, wood- 
land, bare soil, and urban usage. 

There are two types of requirements for land cover classification 
for the RPC 's "280" project. Aggregate land cover statistics for 

water supply, urban storm flow, and estuary pollution are required 
for the three large watershed prototype areas which comprise 19% of the 
total area and which typify critical water quality problems in the area. 
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The statistics derived for each of these relatively large watersheds will 
be used as inputs to runoff models. The estimation of the direct runoff 
is crucial to the determination of total pollutant loads and for the es- 
timation of the first flush phenomena. Since the areas of the prototype 
watersheds are relatively large (>10,000 acres) , the requirement for 
class accuracies of 5% does not impose a severe constraint. Subsequent 
alteration of the signatures has resulted in acceptable accuracy from 
computer classification of LANDSAT for the larger hydrologic basins. 

The second requirement is far more stringent and affects a far 
greater proportion of the area, about 80% of the total. During the 
early part of June, 19.77, it was indicated by RPC that the LANDSAT data_ 
could be given a much more extensive role than orginally planned if the 
accuracies for small areas could approach the requirements indicated in 
Figure 3 . 

High accuracies are required for the small areas since the pol- 
lutant loads are estimated from correlations with in-situ measurements 
related to specially selected small, homogeneous areas. The pollutant 
loads are directly estimated from correlations with the in-situ ground 
sampling sites illustrated in Figure 4. It is estimated that correla- 
tions would exhibit one sigma variance of the order 10-20% given per- 
fect information on land cover. In order to minimize the error con- 
tribution from land classification, this requires that small area land 
classification errors of less than 5-10% be maintained. 

Techniques for improving the small area results were required in 
order to apply LANDSAT data. 



figure 4 
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V. APPLICATIONS OF OPTICAL ANALYSIS TECHNIQUES TO LAND COVER CLASSI- 
FICATIONS FOR NON-POINT SOURCE POLLUTANT LOAD ESTIMATION 

During June, 1977, an Initial presentation was made by ECOsystems 
personnel to the responsible planners and engineers at the RPC. The pre- 
sentation described the utility of LANDS AT data for the determination 
of hydrologically important land cover classes for single season and 
seasonally changing features. Demonstrated accuracies for classifica- 
tion of 95% (5% error) for small watershed areas, typical of those of 
importance to the Baltimore RPC area, were of immediate interest. The 
correlation of the LANDSAT derived information with important hydro- 
logic parameters governing overland flow was also demonstrated. These 
included the determination of land use classifications and the deriva- 
tion of Manning’s "n" from satellite imagery as performed in earlier 
phases of NASS- 30539. 

In order to confirm the utility and accuracy of the optical 
analysis technique for use in non-point sources pollution estimation, 
it was mutually decided to perform a demonstration test on some small 
watersheds located throughout the Baltimore region. The Rhode River 
watershed was selected because of the computer classification data ef- 
fort and ground truth data available. Two other watersheds, in Carroll 
and Howard County respectively, were also selected to prove the extend- 
ability of results. 

5. 1 Demonstration of Land Cover Classification using Optical 
Overlay for Small Test Areas in the Baltimore Region 

The optical overlay procedure employed determines land cover 
classification from LANDSAT using USGS topographic maps as base maps 
and utilizes aerial photographs for ground truth verification. 



II I 


II II I II 


I II I 1 1 ■ 1 1 II 


-*18- 

Three-band combination color diazos are used in the analysis of 
land cover. Spectral Band 5 (0.6-0, 7 )Jm) is encoded in magenta. Band 7 
(0. 8-1.1 ym) is encoded in cyan and Band 4 (0.5-0. 6 ym) in yellow. 

The colors employed are selected to maximize visual discrimination among 
three bands. The cost of a diazo composite, exclusive of the price of 
the original imagery, is about one dollar. 

Since the diazo film is not color stable when exposed to high in- 
tensity projection lights for extended periods, the composites were 
rephotographed on color slide film. The recent introduction of new 
colorfast diazo film, however, can obviate this step in the future. 

Visual interpretation of LANDSAT imagery is- significantly im- 
proved by the use of multi- temporal imagery. The accurate identifi- 
cation of such themes as agricultural usage depends upon the observa- 
tion of their behavior throughout the growing season. LANDSAT land 
cover analysis was performed using multi-temporal images. 

The Ecosystems Universal Image Analyzer, shown in Figure 5, 
was used. Projected images are reflected onto a work table from an 
overhead mirror. Electronic controls at the table allow changes in 
magnification and position of the image. Typically, the projected 
image is overlayed on a map of known scale, such as the USGS 71/2 
minute quad sheet (1:24000 scale). 

The boundary of the watershed of interest is traced from the 
USGS map. Land cover is interpreted from the reflectance charac- 
teristics of the MSS bands, as determined by the color coding of the 
composite. Dense vegetation, for example, gives very low reflectance 
in Band 5. In the projected image, therefore, areas of dense vegeta- 
tion will appear a dark reddish color; this is due to the magenta 
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coding of Band 5 in the diazo composite. Similarly water will appear 
black, urban areas will generally be white to light blue, grass and 
pasture will be a lighter red or orange, bare soil will be light blue 
to blue, and agricultural areas will vary between blue and light red 
depending upon the stage of growth of the crops. When the total 
classification is complete, the areas of each class is planimetered, 

5 . 2 Land Classification of the Rhode River Estuary Test Area 

The optical classification procedure was applied to the Rhode 
River Estuary using LANDSAT imagery from April and October, 1973, and 
February, 1974. The Rhode is a tributary of the Chesapeake Bay lo- 
cated in Anne Arundel County, Maryland as shown in Figure 6. The test 
area covers about 3200 hectares and is predominantly rural with some 
residential development. Other land uses present include forests, 
grass/pasture and agricultural areas. The 1:24000 USGS topographic 
map was used to establish the ridge line of the basin. 

Color infrared aerial photography at 1:60000 scale provided by 
the RPC served as the ground truth. These photos were taken in 
September, 1970, but were deemed accurate for the purposes of this 
analysis as the area was not subjected to any significant developnent . 

Black and white 1975 LANDSAT originals were combined into diazo 
color composites with Band 4 in yellow. Band 5 in magenta, and Band 7 
in cy an . 

Nine small training areas were isolated and an overlay of their 
location was made. The sequential color of each training area was 
noted. The image with the highest discriminability for each class 
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was noted and a preliminary classification of the whole basin was 
made. Successive images were then used to refine the initial clas- 
sification by the multi-temporal overlay procedure. 

The area of each land cover class was then measured by plani^ 
meter and compared to a ground truth map derived from the aerial photo- 
graphs, Table I shows the results of this comparison. The average 
weighted inventory error was 5.5%. 

These results are consistent with earlier findings in that the 
inventory errors are small except in classes with very small absolute 
areas, e.g. , cropped areas, developed areas and surface water. The 
accuracies were commensurate with the RPC requirements and represented 
an improvement over computer-only classifications. 

5 . 3 Land Classification of the Howard County Watershed Test Area 

The Howard County site shown in Figure 7 comprises 4A2 hectares 
of land; the primary land cover class is grass /pasture. About three 
percent of the basin has been developed. It drains into the Little 
Patuxent River, which in turn flows into the Chesapeake Bay. 

The procedure applied in the analysis of this watershed was the 
same as that described for the Rhode River. Because of its small 
size, no sub-site analysis was attempted. 1:60,000 color infrared 
aerial photographs were lised for ground truth. 

The LANDSAT images used were from April and October, 1973, and 
February, 1974. By sequential projection of each scene, it was de- 
termined that only three cover classes were present: forest, pasture/ 
grass, and residential. Forest was segregated using the February 
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TAHLE I 


RHODE RIVER WATERSHED 

RESULTS OF SATELLITE LAND USE CLASSIFICATION 


2 OF TOTAL AREA 2 OF TOTAL 

AREA (MEASURED (HECTARES) AREA (MEASURED 
FROM AERIAL) (AERIAL) FROM LANDSAT) 



AREA 

(HECTARES) 2 

(LAKDSAT) ERROR 


WEIGHED AVERAGE ERROR - 5.5 
































FIGURE 7 


HOWARD COUNTY SITE 
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image. Identification of the other land cover classes was accomplished 
by using the other images in a multi- temp oral sequence. 

Measurements of the land cover classes were performed and coiiqjared 
with the measurements of comparable classes from the aerial photographs. 
The results are shown in Figure 8 and Table II. The weighted aver- 

age error was 2.4% which met the classification criterion established 
by the RPC. 

5 . 4 Land Classification of the Carroll County Test Area 

The Carroll County watershed is located in an agricultural area 
west of Westminster, Maryland. This site, shown in Figure 9, is 
predominately agricultural. Of its 650 acres, about 600 are in fields 
and pastures; the remaining acres are forest. The residential area 
within the watershed is extremely small. The average size of fields 
within the site area is about nine acres. 

The satellite images were prepared as composites with Band 4 in 
yellow. Band 5 in magenta, and Band 7 in cyan. 

The October image was projected and overlayed on a USGS topo- 
graphic map to set the scale. The color of each individual field 
was recorded. 

This procedure was subsequently repeated for the February, 1974 
image. As noted earlier, no ground truth existed for 1973; therefore, 
an absolute comparison of the satellite and aerial images was not pos- 
sible, since fields in the area are rotated from year-to-year. The 
analysis was continued, however, realizing that changes in the field 
content from 1973 to 1975 might be the source of error. 
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FIGURE 8 

HOWARD COUNTY TEST SITE 
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TABLE II 

HOWARD COUNTY SITE 

RESULTS OF SATELLITE LAND USE CLASSIFICATION 


LAND USE CLASS 

% OF TOTAL 
AREA (AERIAL) 

ARE A- AERIAL 
(hectares) 

% OF TOTAL 
AREA (LANDSAT) 

AREA-LANDSAT 

(hectares) 

FORESTS 

35.5 

156.9 

36.1 

159.6 

FIELDS 

61.5 

271.8 

60.3 

266.5 

DEVELOPED 

3.0 

13.3 

3.6 

15.9 

WATER 



- 


- 

- 



WEIGHTED AVERA.GE ERROR = l,k% 
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In examination of the two LANDSAT images, it was noted that two 
distinct classes of fields were present. Roughly half of the fields 
underwent a distinct change in color from October to February, i.e., 
several fields appeared bright orange (indicating dense vegetation) 
in October and light blue (indicating soil or sparse vegetation) in 
February. The remainder of the fields remained constant during the 
season. The former class represented cropped areas while the latter 
were pastures. The results of the classification are given in 
Table III and indicated a weighted error of 5% which met the criterion 
established . 

5 . 5 Conclusions 

From the results achieved for the three test basins, it was ap- 
parent that overall accuracies from optical analysis are satisfactory 
for non-point source pollutant assessment. The average inventory er- 
ror of the three sites was less than five percent. This represents 
a Substantial improvement over the twenty percent errors typical in 
computer-only analysis. These results permit the computation of sur- 
face related hydrologic parameters at desired accuracies. 

The three hydrologic parameters of importance to the pollutant 
load estimates are: surface friction of the overland flow areas of 

the watershed (Manning's “n") , estimation of the seasonal values of 
the ratio of direct-to-delayed runoff due to a given storm event, and 
the impervious /pervious fractions of the watershed. All three of 
these parameters play key roles in the determination of non-point 
source pollution estimates for the non-prototype areas and are amen- 
able to direct estimation from satellite imagery. 
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TABLE III 

CARROLL COUNTY SITE 

RESULTS OF SATELLITE LAND USE CLASSIFICATION 


land COVER 
CLASS 

% OF TOTAL 
AREA 
(aerial) 

AREA 

(hectares) 

AERIAL 

% OF TOTAL 
AREA 

(LANDS AT) 

AREA 

(hectares) 
LANDS AT 

% 

ERROR 

CROPPED 






AREA 

44. 6 

116.0 

42.1 

109.5 

5.6 

PASTURES 

55,4 

144,0 

57.9 

150.5 

4.5 


WEIGHTED AVERAGE ERROR = 5% 
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Optical LANDSAT image analysis was also cost effective. Direct 

comparison of the cost per unit area was not possible since detailed 

records of the alternatives were not available from the RPC. The pro- 

2 

jected costs from the ECOsystems experience of $0. 50/km would yield a 
total cost of $2500 for the RPC area, a figure which was well received 
by the RPC personnel. 

There is an inherent interest at RPC to use LANDSAT type infor- 
mation in a simple and direct way. So far the problems of limited ac- 
curacy, cost of acquisition of ground truth, and difficulty in signa- 
ture extension using conventional LANDSAT computer classification tech- 
niques have limited the utility of the LANDSAT data. The application 
of the simple optical techniques appears to have potential for upgrading 
the utility of LANDSAT. An explanation of how to interpret LANDSAT im- 
ages using these procedures was made to the RPC personnel. The optical 
technique permitted direct interpretation by the multi-disciplinary RPC 
team with enhanced information return. The low cost of the Images per- 
mits expanded analysis of seasonal land cover which closely matches the 
requirement for incorporation into the modeling of pollutant loads as- 
sociated with critical runoff events. 

Based on the results achieved for the test sites and this posi- 
tive reception by the RPC, the complete land cover analysis and pol- 
lutant load assessment of the Magothy River was undertaken. The pro- 
cedure applied and the conclusions reached are reported in the next 


section. 
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VI. ESTIMATION OF NON-POINT POLLUTANT LOADS FOR THE MAGOTHY BASIN 

USING LANDSAT DERIVED DATA 

The estimation of non-point source pollutant loads from represen- 
tative watersheds is a key element in RPC’s approach to areawide waste- 
water management for the Baltimore Region. The anticipated high cost 
of modeling each of the more than 100 water quality districts comprising 
the region required that the RPC select a limited number of "prototype" 
basins which were adjudged to be representative of the region and from 
which the total non-point source loads of the non-prototype areas could 
be estimated. 

One of the three prototype regions selected by the RPC was the 
Magothy River basin located in close proximity to the ECOsystems 
facility at Gambrills, Maryland. The Magothy is a 8600 hectare 
watershed which drains into the tidally dominated Magothy River as 
shown in Figure 10. 

The ECOsystems model ECHOS* was used to calculate annual pollutant 
loads for the Magothy, The model computes overland flow from indivi- 
dual rain events and stores the output as flow versus time. The pro- 
cess is repeated for each sub-basin. Subsequently, the model routes 
the output of the appropriate sub-basin or sub-basins through suc- 
cessive channels until the basin outlet has been reached. Tabular re- 
sults are prepared giving the time since the start of rainfall, rain 
and discharge rates, and cumulative rain and discharge. Total dis- 
charge for the month is computed and converted to pollutant loads using 
loading factors determined for the existing land use categories. 


*See Appendix A for a description of the model. 
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table IV 

COLOR COMPOSITE ANALYSIS; MACOTHY RIVER BASIN 
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Table IV indicates that no single image shows agricultural areas 
in marked contrast with the other classes. It was determined by examina- 
tion of all the available scences that agricultural areas followed a 
unique progression of colors: dark pink in April, light blue to white 
in October, and orange/light blue in February. All areas exhibiting 
this behavior were then classified as agricultural. 

Residential areas were not strongly separable from pasture/grass 
in any scene. This is due to the low density (2-4 units per acre) of 
most residential development present. Some of the more densely de- 
veloped residential areas could be consistently identified; most ap- 
peared to be the same as the pasture/grass areas in all scences. The 
residential class was determined from aerial photography. The use of 
a hybrid satellite/aerial remote sensing system to derive pollution 
loads represents a necessary compromise. 

On the basis of this technique, the land classification map shown 
in Figure 12 was prepared. The entire Magothy River watershed was also 
mapped and classified from the aerial photo. Because of-, the age of the 
aerials (September, 1970), imagery from 1975 was acquired from the Baltimore 
RPC and the ground truth map updated. Alterations were required in 
ten areas reflecting the post 1970 residential and/or commercial develop- 
ment . 

The areas of each land cover class were computed from the aerial and 
hybrid maps for each of the seventeen subbasins. The results are re- 
ported in Table V. The overall classification accuracy of 2.3% is con- 
sistent with the other watersheds analyzed and well within the require- 
ment of the RPC. The land use classifications was also used to develop 
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table V 


LANDSAT a AERIAL CLASSIRCATION 

AKEA, hectares 


OF THE MAGDTHY RIVER 


SUBBASIN 

FOREST 

RESIDENTIAL 

COMMERCIAL, 

INDUSTRIAL 

WATER 

EXTRACTIVE 

AGRICULTURE 

PASTURE, 

CRASS 

TOTAL 

1 


25.6 

- 

18.6 

- 

- 

- 

116.2 

2 


255.6 

- 

- 

- 

- 

74.4 

669.2 

3 


204.5 

- 

- 

- 

- 

- 

485.6 

4 

237.0 

25.6 

- 

- 

- 

- 

- 

262.6 

5 

474.0 

227.7 

- 

- 

- 

- 

55.7 

757.4 

6 

1279.2 

697.1 

69.7 

5.8 

41.8 

190.5 

86.0 

2370.1 

7 

KM 

51.1 

- 

- 

- 

- 

- 

2 32,3 

B 


209.1 

13.9 

- 

- 

- 

39.5 

697,0 

9 

106.8 

160,3 

34.9 

- 

- 

- 

- 

302.0 

10 

97,6 

41.8 

11.6 

- 

- 

- 

- 

151.0 

11 

260.1 

48.9 

20.9 

- 

- 

- 

65.1 

395.0 

12 

188.2 

165.0 

16.3 

- 

- 

- 

95.3 

464.8 

13 

300.8 

178.9 

10.5 

- 

- 

- 

67,4 

557.6 

14 

260.3 

264.8 

25.6 

- 

- 

- 

53,4 

604.1 

IS 

241.7 

162.6 

2.3 

- 

- 

23.2 

104.6 

534.4 

TOTAL 

4753.7 

2718,6 

205.7 

24.4 

41.8 

213.7 

641.4 

8599.3 

ERROR 

-1.51 

- 

-1.9Z 

-9.6X 

-5.2Z 

-8.9Z 

+18. 9Z 



WEIGHTED AVERAGE ERROR - 2.3Z 
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the required hydrologic inputs to the ECHOS model and to assess the 
pollutant load conversion coefficients for the 17 s ub -water she ds . 

6 . 2 Determination of the Rainfall and Runoff Inputs to the Model 

Storm generated, non-point source load estimation utilizes the 
correlation between land cover generated pollutant loads and the volume 
of runoff. Sediments and sediment-entrained pollutants are carried 
overland by stormwater runoff and deposited into the stream channels. 
After deposition in the channels, they are transported and altered by 
channel flow processes to the stream outlet. The process of pollutant 
delivery to the channel outlet can be simulated by computer models with 
routines for overland and channel flow, such as the ECHOS model developed 
under NAS 8- 305 39. 

The ECHOS model requires calibration with a small number of actual 
rainfall and runoff records. Though the Magothy River is not perma- 
nently gaged, a limited amount of on-site data was collected by the RPC. 
Historical and current rainfall records were also available for nearby 
weather stations. 

On-site data were taken for various months during 1977 at three 
sites: Bay Hills, Sevema Park, and Cape St. Clair. Measurements 

were taken for a total of 80 rainfall events. However, sufficient 
runoff records existed for only four events at Cape St. Clair; 24 
February, 18 March, 2 April, and 25 July, 1977 from which to con- 
struct an accurate hydrograph. Rainfall hyetographs and runoff hy— 
drographs for these dates are graphed in Figures 13a, b, c, and d. 

These graphs show the actual runoff response of Magothy sub-water- 
sheds to actual rain events. 



RAIN RA1ES 


FIGURE 13a 





700 aoo 900 1000 1100 1200 1X0 1400 1500 1600 1X0 

TIfC (HOURS) 


RDHorr 




<■/«( X (10"®) RAIN RATE («/b) 


-40- 

PIC-URE 13c 


rainfall and runoff profiles at cape ST. CLAIR 
DATE: 2 APRIL 1977 
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Monthly averages of number of storm events, volume, duration, inten- 
sity and time between storms were computed from an analysis of twenty- 
five years of hourly rainfall records from the nearby NOAA weather 
station at Baltimore Friendship Airport. Typical rainfall data avail- 
able from NOAA is shown in Table VI. The raw data were processed by 
computer to produce a statistical summary. The monthly summary pre- 
sented in Table VII gives the intensity, volume, duration and time be- 
tween storms and includes the average, minimum, maximum, standard de- 
viation, variance and C.V. for each parameter. Table VIII summarizes 
the data from the Table VII and lists the selected parameters of the monthly 
average rain events to be used in the simulation. 

Direct runoff from a storm of a given volume, magnitude and* aver- 
age intensity depends directly upon the distribution of the rain rates 
over time. Discharge rates and volumes vary with the "shape” of the 
rain input. Average peak rain rate and average time to peak rate were 
derived for each month using a large, fifty percent sample of the storms 
from 1948-1975 and recording the peak rate and time to peak rate. Ihe 
sample comprised from 100 to 150 storms, depending on the month for the 
period of record, and were representative of the actual rainfall behavior. 

The temporal rainfall profiles were approximated by a series of 
discrete intervals with total volume, intensity, peak rate, and time to 
peak equal to the average values adjusted to match the durations and 
volumes specified in Table VII. The June 1975 profile is shown in 
Figure 14 along with its discrete approximation as an example of this 
procedure. The average number of storms expected during each month 
was computed from the average time between storms and is also given in 


Table VIII. 
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table VI 

HOURLY RAINFALL DATA; BALTIMORE AIRPORT GAGE FOR THE 
PERIOD JANUARY - JUNE 1975 
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9ft 

ft 

0 

ft 

ft 

ft 

ft 

0 

ft 

ft 

0 


1404 70 

2/17/79 

0 

0 

0 

0 

0 

10 

9 

0 

1 

1 


ft 

0 

0 

ft 

ft 

ft 

ft 

0 

0 

0 

0 

ft 


140470 

2/14/79 

0 

0 

1 


0 

0 

1 

0 

0 

0 


0 

0 

ft 

0 

ft 

0 

ft 

0 

0 

0 

0 

0 


1104 70 

2/23/79 

0 

0 

0 

0 

0 

0 

0 

2 

27 

27 

10 

0 

0 

ft 

0 

0 

ft 

ft 

ft 

0 

ft 

ft 

ft 


1404 70 

2/24/79 

2 

1 

t 

0 

0 

0 

0 

10 

3 

3 


0 

1 

3 

9 

2 

1 

1 

1. 

3 

0 

ft 

ft 


140470 

3/ 1/79 

0 

0 

0 

0 

0 

0 

• 

ft 

ft 

0 


0 

0 

0 

ft 

ft 

0 

ft 

0 

0 

0 

ft 

ft 


140470 

3/ 7/79 

0 

0 

t 

0 

0 

0 

0 

0 

ft 

ft 


0 

1 

3 

ft 

0 

0 

ft 

ft 

ft 

ft 

ft 

ft 


140470 

3/10/79 

0 

0 

0 

• 

0 

# 

0 

0 

ft 

ft 


ft 

2 

2 

4 

0 

0 

0 

0 

0 

ft 

0 

ft 


140470 

3/12/79 

0 

0 

0 

0 

e 

0 

ft 

ft 

2 

ft 


Ift 

ft 

1 

14 

L 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


1404 70 

3/13/79 

0 

2 

i 

0 

0 

0 

0 

ft 

0 

ft 


ft 

1 

ft 

0 

ft 

ft 

ft 

0 

ft 

0 

0 

ft 


140470 

3/14/79 

• 

14 

9 

4 

0 

4 

1 

ft 

ft 

1 

24 

22 

9 

12 

3 

ft 

0 

ft 

0 

ft 

ft 

3 

ft 


110470 

3/19/79 

0 

0 

0 

0 

3 

19 

1ft 

T 

24 

9 

IT 

49 

39 

44 

4 

0 

0 

0 

0 

0 

0 

ft 

ft 


1404 70 

1/24/79 

0 

0 

0 

0 

4 

f 

4 

31 

7 

0 


0 

0 

0 

9 

19 

12 

0 

ft 

ft 

ft 

ft 

ft 


1404 70 

3/29/79 

1 

0 

0 

0 

0 

3 

0 

ft 

1 

4 


0 

ft 

0 

ft 

ft 

0 

ft 

ft 

0 

ft 

0 

ft 


1404 70 

3/30/79 

11 

2 

10 

14 

0 

4 

0 

0 

ft 

ft 


0 

ft 

0 

ft 

ft 

0 

0 

0 

ft 

ft 

0 

ft 


1404 70 

4/ 1/79 

0 

0 

0 

0 

0 

0 

0 

0 

ft 

ft 


0 

0 

ft 

0 

ft 

0 

ft 

0 

ft 

ft 

ft 

ft 


1404 70 

4/ 3/79 

0 

0 

0 

0 

1 

2 

1 

2 

24 

0 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

0 

ft 

ft 

0 

ft 


140470 

4/19/79 

0 

4 

• 

9 

3 

9 

1 

1 

1 

4 


0 

ft 

9 

2 

0 

ft 

ft 

0 

ft 

ft 

ft 

ft 


1404 70 

4/14/79 

0 

0 

0 

0 

• 

0 

ft 

ft 

0 

0 


2 

ft 

ft 

0 

0 

0 

ft 

0 

0 

ft 

0 

ft 


140470 

4/19/79 

0 

0 

0 

0 

0 

0 

ft 

ft 

ft 

0 


19 

2 

ft 

ft 

ft 

0 

0 

0 

0 

ft 

ft 

ft 


1404 70 

4/23/79 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ft 


ft 

ft 

ft 

0 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


1404 70 

4/24/79 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

ft 

0 

9 

ft 

1 

10 

12 

11 

1 

0 


140470 

4/29/79 

4 

2 

0 

Q 

D 

0 

0 

ft 

0 

• 


ft 

ft 

4 

1 7 

10 

• 

49 

24 

13 

13 

4 

0 


140470 

4/24/79 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


ft 

ft 

ft 

ft 

0 

ft 

ft 

0 

0 

0 

1 

2 


1404 70 

4/29/79 

0 

0 

0 

12 

3 

1 

2 

0 

ft 

0 


ft 

ft 

0 

0 

9 

ft 

0 

0 

ft 

D 

0 

ft 


140470 

9/ 1/79 

1 

1 

4 

14 

17 

IV 

4 

2 

1 

ft 
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ft 

ft 

0 

ft 

0 

0 

ft 

0 

ft 

ft 

ft 


140470 

9/ )/79 

0 

0 

0 

0 

0 

0 

ft 

D 

ft 

0 


D 

0 

ft 

ft 

ft 

0 

ft 

0 

ft 

ft 

3 

4 

19 

140470 

9/ 4/79 

14 

14 

21 

12 

u 

33 

11 

7 

9 

0 


0 

0 

ft 

ft 

22 

i 

0 

ft 

9 

27 

0 

ft 


140470 

3/ 5/79 

0 

1 

0 

0 

0 

0 

0 

0 

ft 

0 


ft 

0 

ft 

0 

ft 

0 

ft 

0 

ft 

0 

0 

ft 


140470 

9/ 4/79 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

ft 

0 

ft 

1 

47 

0 

0 

ft 

ft 

0 

ft 


1474 70 

9/12/79 

0 

0 

0 

0 

0 

g 

0 

ft 

0 

0 


0 

ft 

ft 

9 

ft 

40 

1 

ft 

0 

0 

tl 

14 


140470 

9/11/79 

0 

0 

• 

0 

0 

• 

0 

0 

ft 

0 


ft 

ft 

0 

3 

0 

ft 

0 

2 

0 

0 

0 

ft 


140470 

9/11/79 

0 

0 

0 

0 

0 

0 

0 

ft 

ft 

0 


ft 

ft 

ft 

0 

ft 

0 

0 

0 

ft 

ft 

1 

4 


I4M70 

9/14/79 

11 

3 

1 

1 

1 

1 

2 

1 

0 

0 


0 

0 

9 

0 

0 

ft 

ft 

ft 

ft 


0 

0 


1404 70 

9/22/79 

0 

0 

■ 

0 

0 

0 

ft 

ft 

0 

ft 


0 
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ft 

0 

ft 

0 

39 

21 

ft 

0 

0 

0 


140470 

9/24/ 79 

0 

0 

0 

0 

0 

0 

0 

ft 

ft 

ft 


112 

4i 

0 

ft 

0 

0 

ft 

0 

0 

ft 

0 

D 


14^70 

9/10/79 

0 

Q 

0 

0 

0 

0 

ft 

ft 

ft 

0 


1 

14 

ft 

ft 

ft 

0 

0 

0 

0 

9 

ft 

0 


140470 

9/JI/79 

0 

0 

0 

0 

0 

0 

ft 

0 

0 

0 


9 

0 

ft 

ft 

0 

0 

ft 

ft 

1 

4 

34 

19 

27 

140470 

4/ 1/79 

1 

0 

0 

0 

I 

4 

10 

9 

ft 

ft 


0 

9 

0 

22 

2 

ft 

0 

ft 

ft 

0 

0 

ft 




43- 


T/^BLE VII 



MONTHLY STATISTICAL SUMMARY, 1948-1975 FOR THE 
BALTIMORE AIRPORT RAIN GAGE. 


KAINFALL statistics ST MOM1HIFOS «IOD Of RICOAOl 

month MUHBEK TOTAt H|N|HUH MAXIMUM AVERA6C STD DEV 

1 

DURATION 2A2. D«ltS100E»04 O.IOOOOOE*01 0«A4000QE»02 O.T44f?AE«Dl O.AVAA3VE«OI 

IHTENSIIT 242. 0«A2tSSTE»0l D«42t*O0E-O2 O«lTS0Olf»OD 0.S1V4ITE-01 0.2S1I4TE-0I 

VOLUME 242« 0.7464AbE»02 O.IOOOOOE’OI 0«ZS700Df «^Ot 0*3aT2T2E»Oti 0.4ll5bOEAOO 

delta 242« 0«202b4SE^Q5 0*4SOOOOE«01 0«740000C*0) 0*S1737AE«02 0*T0VD82E^02 

2 

OLWATIOM 221* O.lAltOOE^DA 0.10000QE»0I 0.440000E»02 0«7V42VBE»Ot 0.tA43tlE«01 

INTENSITV 22B« 0.V7IS54€^01 O.SOOlOOE-02 0«3 BOODSF ♦00 0»42V|VO€-Ot 0*4V701IE-0i 

VOLUME 22A, 0.f767B7E*02 0« 1000001*01 O.IVAOOOE^OI 0«3A4394E^OO 0*434T40E^00 

delta 22B« 0*ll0290E^aS O.SOOOOOE»Oi 0«493))OQE^03 0«79J744E^02 0*T7003lE^02 

3 

OLMATIOfi 24M, 0,iV4400E»04 0.100000E*01 0.350000E»02 0«722477E*01 0.443R28E*0& 

iNTENSITT 24V. 0.I333>4E»02 O.4OOI0OE*02 O.LIOOUJE^OI 0.49S740E*0I 0.7VIK42E-01 

VOLUME 24 V. 0.10444ae*03 0.1O0000C*0l O.424300E^0l 0.3a7027E*00 O.337747E^00 

DELTA 24V. 0.2047 IOE*03 0.400000E*01 0«413VOOE^03 0.T4t4426»02 0.T31712E^02 

4 

OUMATim 302. 0.L7|400E«O4 O.IOQOOOE^QI 0«IIDOOOE»03 0.347530E40I 0.7l8l44E*0l 

INTENSITT 302. 0.143814E902 O«571S00E*02 O.320OO9E^O0 0.4i2l33E*Ol 0«42VSI4E*0l 

VOLUME 302. Q.V324V2E^02 0.1000006*01 O.SVAOOOE^Ol 0.30074VE»09 0«S0i733E400 

DELTA 302. 0.1V3I20E«09 0.4QOOOOE401 0.30S300E403 0.43V470E»02 0.440B31E«02 

9 

OUAATICM 312. 0.141200E«04 0.100000E»01 0.370000E«02 0.91444TE401 O.S54470E«Oi 

iNTtNSlTT 312. 0.217V14E»02 0.9001002*02 0.473002E400 0.4V«490E*0I O.VOISSVE-Ol 

VOLUME 312. 0.10790Sf^03 0.100000E*01 0*3V1000E«OI 0.344977E»00 0.4V3f9VE«00 

DELTA JI2. 0.2O93IOE«O9 0.0 0.394000E»01 0.498049E«02 0.424479E«Q2 

4 

DURATION 2Vl. 0.1U300E404 O.lOOOOOEtOl 0.240000E«02 0.3VV4S4E*01 0.403447E401 

INTENSITY 2V1. 0.31I471E402 0.5001D0E-Q2 0.220000E«01 0. |07039E^00 0.1*47V7E^00 

volume 2Vl. 0.10V41AE^03 0.|iM0O0E-01 0.440300E«01 0«S7400AE»00 O.9970l4E^00 

DELTA 2V1. 0«204r79E»09 O.AOOOPOE^Ol 0.97|>OOE^03 0.710947E*02 O.I3TI34E«^02 

T 

DURATION 294. 0.tt7000E^03 O.IOOOOOE«01 0.343000E»02 0.34V213E401 0.379930E401 

INTENSITT 294. 0.12404»€^02 0.9O0i00E*O2 O.IWQOOE^OI 0.1279TAE»00 0.2002I3E»00 

VOLUME 294. 0.10744lE^03 O.IOOOOOE*Q1 0.94IOOOE»Ol 0.423029E#00 0.7I2471E»00 

DELTA 294. 0.1TV020£^09 3.430000E401 0.3A40006403 0.7I3943E#02 0.742t48E»02 

• 

DURATION 244. O.iOVSOOE^OR O.lOOOOOEtOl 0.3100006*02 0.41L494E*0l 0.4SI2V4E*Oi 

INTENSITT 244. 0.30At04E*02 0.900lOOE*02 0.lOl0O9E*OI 0.1lS092E*00 0.134I9VE*00 

VOLUME 244. O.I2VltOE*03 0* IOOOOOE*OI O.7VS0O0E*01 0.4I34T1E*00 0.093424E*00 

DELTA 244. 0.1Vt009E*09 0.400000€*01 0.4I09U0E*03 0.744340€*Q2 0.74A2V0E*02 

0 

DURATION 224. O.I14700E*04 0,1000006*01 O.43ODOOE*02 0.900RT3E*Ol 0.9T2032C*Ol 

INTENSITT 22V. 0.1TI649E*02 0.0001006-02 0.4071436*00 0.7I010VE-0I 0.1011326*00 

VOLUME 22V. O.V44204E*O2 O.IOOOOOE-01 0.4293001*21 0.4210»4E*00 0.72lS34E»00 

DELTA 22V. 0.2I3V09C*09 0.400000E*OI 0.4340006*03 O.V34003E*02 0.1103ViE*03 

10 

DURATION lOV. 0.tl4900E*04 Q»tOOQOOE*01 0.240aoaC*02 0.4144026*01 Q.944498E«01 

INTENSHT lOV. 0.1002496*02 0.444400E-02 0.30444IE*00 0.9I2V34E-0t 0.9494406*01 

VOLUME IRV. 0.424MVE*02 0.1000006*01 0.49iOOOE*01 0.43444VE*00 0.493903£»00 

delta liv. 0.1V1430E*09 0.4900006*01 O.9O*000E*03 0.1090996*03 O.I041V3E*03 

11 

DURATION 214. 0.1474006*04 0.1000001*01 O.9V0O0OE*O2 0.44VT20E*0l 0.474I34€*0I 

INTENSITY 214. 0.1071316*02 0.9001096*02 0.2990006*00 0.9004746*01 0.4120416*01 

VCLUMf 214. 0.9tV34TE*02 0.1000006*01 0.9140006*01 0.42V420E*00 0.4472196*00 

DELIA 214. 0.201V106*09 0.4900006*01 0.4V9000E*QJ Q.V749426*01 0«10043TE*03 

12 


variance 

0.4t902TE»02 

0.7V44V3E*03 

D.14V3426*09 

0.1244906*04 

0.41710fE*02 
0.247020E*02 
0.204V roEtoa 
0.9V2V476*04 

0.443327E*02 

0.424TWE-02 

0.21V172E*00 

0.93940ZE*04 

0.421Z03E*OZ 

0. 1144846-OZ 
0.Z91T01E*03 
0.4104196*04 

0.30911ZE*02 

0«i29ll9E-02 

0.Z43VV9E*00 

0.3VZ470E*04 

0.1419316*02 
0.3419326*01 
0. 3102446*00 
0.7001276*04 

0. 1413241*02 
0.4001936*01 
0.907419E*0a 
0.9119476*04 

0. 14b0L9E*02 
0.2444A2E*01 
0.7214 746*00 
0.959V31E«04 

0.3909246*02 
0. 1022746*01 
0.9312036*00 
0. 1211716*09 

0.1114126*01 

0.3199496*02 

0.4270446*00 

0.1170976*09 

0.4947306*02 
0.2323126*02 
0.4492496*00 
0. 1001746*09 


COEF-VAR 

0.9109296*00 

0.1302746*00 

0,1330496*01 

0.1019436*01 

0.4137776*00 
0. 1191026*01 
0.1113036*01 
0.9738036*00 

0.9213346*00 

0.199*476*01 

0.1312296*01 

0.9919306*00 

0.1318716*01 

0.1199746*00 

0.1424146*01 

0.1002136*01 

0.1077456*01 

0.1300946*01 

0.1433921*01 

0.9920246*00 

0.1009916*01 
0.1749206*01 
0.1481396*01 
0. 1171196*01 

0. 1074916*01 
0. 1949 346* 01 
0. 1484236*01 
0.9734136*00 

0.1047716*01 

0.1344146*01 

0.1791426*01 

0.1009296*01 

0.1112046*01 
0.1294316*01 
0. 1730196*01 
0.1111896*01 

0.9197306*00 

0.9172496*00 

0.1497326*01 

0.1029416*01 

0.9774946*00 

0.9432016*00 

0.1993206*01 

0.1021906*01 


nuRlTlON 

intlnsiiv 

VCLUM6 


111. 0.1939006*04 0.1000006*01 0.3310006*02 0.1174196*01 0*4424476*01 0.4391996*02 0.7449946*00 

214. 0.9490446*01 0.9001006-02 0.184447E*00 0.442490E-OI 0.1988996-01 0.121777C-02 0.4104216*00 

214. U.I091996*03 O.IOOUOOE-OI O«3990OU6*0a 0.4429436*00 0,9422396*01 0.3141136*00 0.1149116*01 
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TABLE VIII 

SUMMARy OF THE PARAMETERS OF THE AVERAGE MONTHLY RAINFALL FOR 
BALTIMORE AIRPORT GAGE 1948-1975 


1 

! 

DURATION 

hr. 

VOLUME 

cm 

1 

INTENSITY 
cm/hr. ' 

P 

c 

cm/hr. 

■ 

AVERAGE 

NUMBER 

STORMS 

January 

9.12 

7.9 

o 

1.7 

3.37 ; 

8.14 

February 

8.96 

9.8 


2.4 

3.89 

7.63 

March 

i 

7.55 1 

9.9 


2.3 

3.31 

8.78 

April 

6.40 

7.8 


2.6 

2.97 

10.24 

May 

4.93 

8.8 

1.8 

3.6 

2.74 

10.52 

June 

1 

3.52 

9.5 

2.7 

4.8 

1.71 

9.65 

July 

3.32 

0.8 

3.2 

6.1 

1.79 

9.11 

August 

4.18 

2.3 

2.9 

5.8 

1.98 

9.48 

September 

5.39 

10.7 

2.0 

3.6 

2.63 

7.29 

October 

7.62 

11.1 

1.5 

3.2 

3.20 

6.60 

November 

8.58 

10.9 

1.3 

3.1 

3.60 

6.78 

December 

11.12 

12.3 

1.1 

3.8 

3.85 

7.46 1 








RAINFALL RATE (m/hr,) 
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For each month, the model computes the runoff for the average 
storm for that month. The pollutant load is computed for that storm. 
Total monthly pollutant loads are computed by multiplying the load for 
the average storm by the average number of storms for the month. 

While the ECOsystem ECHOS model was designed to utilize the time 
history of runoff, the anticipated RPC data relating pollutant load to 
time history of discharge was not available. The RPC indicated that the 
relationships it had derived for pollutant load to total runoff would 
provide the information desired on total annual pollutant loads. Sub- 
sequent improvements with the time history of runoff should improve 
upon the results. Appendix A gives a detailed listing of the ECHOS 
model with examples of the’ -data input required, 

6 . 3 Calibration of the Runoff Portion of the ECHOS Model 

Calibration of the ECHOS model requires the measurement, deriva- 
tion, or estimation of the hydrologic parameters which permit it to 
most accurately reflect the actual behavior of the basin. Most of the 
parameters required are relatively easy to determine. Average slope, 
area, flow lengths, and channel dimensions are readily available from 
sources such as USGS topographic maps. Other parameters, however, are 
not so well understood or documented. It is particularly difficult, 
for example, to have accurate knowledge of movement of sub-surface 
water or soil moisture content just prior to a rainfall event. 

These parameters were estimated for the Magothy River Watershed 
using the rainfall and runoff records described in the previous sec- 
tion, The calibration was made for the Cape St, Clair sub-basin and 
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the results extrapolated to the remainder of the watershed. A map of 
Cape St. Clair is shown in Figure 15. 

Hanning's "n" was estimated directly from this LANDSAT classifi-*- 
cation by using the relationships given in Table IX. Preliminary es- 
timates of soil moisture parameters were made from the SCS soil survey 
for Anne Arundel County, Maryland and used the SCS land cover rela- 
tionships. These included final (saturated) infiltration rate, soil 
moisture capacity and percolation rate. Antecedent soil moisture was 
estimated using the following relationship adapted from the US Weather 
Bureau: 


where: 


ASM 


(n) 


= P 


(n) 


+ .9ASM 


(n-1) 


ASM^ = Antecedent moisture on day (n) 

= Precipitation on day (n) 

ASM, = Antecedent moisture on day (n-1) 

(n-i) 


The initial value of the antecedent soil moisture (ASM) was as- 
sumed to be zero immediately after the longest dry period occurring 
from ten to twenty days prior to the day of the test event. Values 
of ASM were then computed for each day , up to the day of the test 
event. The final value was used as an initial estimator for the ASM 
for each event. 

Once the initial value of all parameters required by the model 
was determined, a series of calibration runs were made for the March, 
April, and July events. Discrepancies in the runoff rate recorded for 
the February event precluded its use for calibration. 




SOURCE: BALTIMORE REGIONAL PLANNING COUNCIL 
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TABLE IX 

MANNING’S ROUGHNESS COEFFICIENT FOR OVERLAND 
FLOW FOR VARIOUS SURFACE TYPES 


WATERSHED SURFACE 

MANNING’S "N" 

Smooth Asphalt 

0.013 

Concrete (Trowel Finish) 

0.013 

Rough Asphalt 

0.016 

Concrete (Unfinished) 

0.017 

Smooth Earth (Bare) 

0.018 

Firm Gravel 

0.020 

Cemented Rubble Masonry 

0.025 

Pasture (Short Grass) 

0.030 

Pasture (High Grass) 

0.035 ' 

1 

Cultivated Area (Row Crops) 

0.035 

Cultivated Area (Field Crops) 

0.040 

Scattered Brush, Heavy Weeds 

0.045 

Light Brush and Trees (Winter) 

0.050 

Light Brush and Trees (Summer) 

0.060 

Dense Brush (Winter) 

0.070 

Dense Brush (Summer) 

0.100 

Heavy Timber 

0.100 

Idle Land 

0.030 

Grass Land 

0.032 
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Results of the first runs showed that the initial percolation 
rates chosen were too large, causing the upper soil layer to drain too 
quickly. This parameter was adjusted and subsequent rims made until a 
satisfactory fit was realized. An example of the goodness-of-fit of the 
technique is shown in Figure 16 for the 25 July 1977 event. 

At this stage in the calibration, all but one of the parameters 
were fixed. Final calibrations were made by adjustment of antecedent 
soil moisture. It was anticipated that some refinement might be neces- 
sary to compensate for the fact that average, -rather than peak rainfall data 
were to be used as inputs to the model. The average rain intensity 
derived from the statistics of the Airport gage were' generally lower 
than the three test events at Cape St. Clair. 

In order to correctly adjust the antecedent soil moisture, it was 
necessary to compute the fraction of runoff expected from rainfall for 
each month. Values for the months of the test events could be calcu- 
lated directly as shown in Table X. Alternately and for the general 
case of limited data, the SCS curve numbers procedure using LANDSAT 
derived land cover was used to test the validity for ungaged basins. 

For the remaining months, the average value of 1,78 was used. This 
fraction was applied to the average rainfall volumes by month to com- 
pute expected monthly runoff. Figure 17 shows the computed monthly 
runoff volume as a % of the rainfall for a typical sub-basin. Dividing 
Creek, 

The final calibration runs were made by adjusting antecedent 
moisture. Dividing Creek was used due to its homogeneous composi- 
tion, relatively large size, and monostream configuration. 



FIGURE 16 

RESULTS OF THE CALIBRATION RUN FOR CAPE ST. CLAIR WATERSHED FOR THE 

25 JULY 77 EVENT 

watershed: capstclr 

EVENT OF! 25JUL77 

♦ INDICATES OBSERVED DISCHARGE 

* INDICATES SIMULATED DISCHARGE 
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TABLE X 


RUNOFF AS A % OF RAINFALL FOR THE CAPE ST. CLAIR 
SUBBASIN CALIBRATION EVENTS 



RAINFALL 

VOLUME 

(cubic meters) 

DISCHARGE 

VOLUME 

(cubic meters) 

RUNOFF AS 
A % OF 
RAINFALL 

FEB 

1982 

20.0 

1.0% 

MAR 

1353.8 

32.65 

2.0% 

APR 

5560 

72.3 

1.3% 

JULY 

3771.2 

87.3 

2.8% 






AVERAGE PERCENT DISCHARGE = 1.78% 
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Three rxans were inacle foir each mont^ at different va?uues of antecedent 
moisture. The volume of runoff was recorded and plotted on a graph of 
discharge versus antecedent moisture. Interpolated values of ASM at 
zero error were then read directly from t:he graph. The results of re- 
peating this process for each month are given in Figure 18'. This 
curve is indicative of rational hydrologic behavior. Lower values of 
soil moisture are evidenced in months with higher evapo transpiration 
and solar radiation and vice-versa. The antecedent moisture figures 
derived for Dividing Creek were assumed to hold for all other sub- 
basins . 


6 . 4 Derivation of Pollution Loading Factors for Magothy River 
Basin 

Non-point source pollutant estimation from land cover classifica- 
tion involves the use of empirical loading factors which relate pol- 
lutant loads to the land use classes, generally: 

m n 

PL. = I Z (LU. X LF. .) 

^ j=i i=i ^ 

where ; 

PL. = Pollutant load, volume/time 

LU^ = Area of land use type i, area 

LF^. = Loading factor for land use type i and pol- 
lutant type j, volume/unit area/ time 

The available data supplied by the RPC, indicated that the estimation 
of pollutant load from land cover measured by LANDSAT was the prefer- 
able technique. Pollution assessment from streamflow, precipitation. 
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and sediment load are generally more applicable to very large areas 
and do not adequately reflect the response characteristics of small 
basins. The general form of the sediment-pollutant coefficients 

would not account for loading differences among the smaller Magothy 
sub-basins. The use of land class - pollutant loading factors directly 
uses the surficial and temporal characteristic of the watershed and is 
therefore directly amenable to remote sensing analysis. 

The applicability of loading factors derived for one geographic 
area to another locale is not certain. It is far preferrable to use 
loading factors developed in the region under study. The availability 
of extensive water quality sampling at the RPC provided ECOsystems with the 
necessary data. 

With the assistance and at the direction of the RPC, two such sets of 
loading coefficients were investigated. One was developed for water- 
sheds in North Carolina and the other by the Washington Council of 
Governments for Maryland and Virginia basins. The data for the pol- 
lutants of interest the RPC are shown in Tables XI and XIT . Ti'.ese data 
sets were judged suitable for the Magothy by the RPC because of their -geo- 
graphic proximity and because they showed a reasonable numerical con- 
sistency for the individual pollutants. 

The loading factors were stored on computer disk files to be read 
into the ECHOS model. After the cumulative monthly discharge was com- 
puted for a sub-basin, the program calculates a weighted average pol- 
lutant concentration based upon the fraction of the sub-basin in each 
land use class. These factors are then multiplied by the discharge to 
yield total pollutant loading. This? procedure is repeated for each 



POLLUTANT LOADING 


LAND USE 

SUSPENDED 

SOLIDS 

(mg/1) 

BOD 

(mg/1) 

COD 

(mg/1) 

Shopping Ctr. 

74 

No Data 

54.66 

Strip 

Commercial 

98 


111.53 

High Rise 
Residential 

34 


96.28 

Townhouse/ 
Garden Apts. 

115 


73.46 

Med. Density 
Residential 

42 


95.88 

Large Lot 
Residential 

142 


115.15 

Construction 

763 


89.35 

Forest 

No Data 


72.40 

Conventional 
Till Agri. 

1681 


853.4 

Minimum 
Till Agri. 

133 


86.24 

Pasture 

544 


549.00 


-Liia.jjtii AjL 


FACTORS - WASH/ COG 


1 TOTAL 
KJELDAHL 
NITROGEN 
(mg/1) 

TOTAL 

NITROGEN 

(mg/1) 

TOTAL 

PHOSPHOROUS 

(mg/1) 

1.81 

2.43 

0.20 

1.33 

1.76 

0.27 

0.87 

1.20 

0.2'2 

1.12 

1.49 

0.29 

1.94 

2.42 

0.38 

2.76 

3.64 

0.35 

1.61 

1.76 

0.26 

i.28 

1.30 

0.12 

8.90 

13.42 

3.18 

1.93 

3.52 

1.14 

4.97 

5.16 

1.06 


I 

I 




















TABLE XII 


POLLUTANT LOADING FACTORS - RPC - TRIANGLE 5 


LAND USE 

% IMPERVIOUS 
AREA 

SUSPENDED 

SOLIDS 

CONCENTRATION 

(mg/l)_ 

BOD 

CONCENTFATION 
(ms/ 11 

COD 

CONCENTRATION 

(mg/1) 

TOTAL KJELDAHL 
NITROGEN 
CONCENTRATION 
... (mg/1) 

TOTAL PHOSPHOROUS 
CONCENTRATION 
(mg/1) 

Rural 

4 

167 

7,0 

38.2 

0.71 

0.41 

LAC 

12 

239 

7.6 

61,1 

1.03 

0.47 

LAR 

16 

293 

7.7 

76.1 

2.32 

0.49 

HAR 

32 

751 

12.2 

140.6 

2.34 

0.68 

HAC 

35 

633 

13.2 

172.8 

2.22 

0.60 

CBD 

80 

179 

15,3 

178.2 

2.97 

0.46 


LAC = Low Activity Commercial 
LAR = Low Activity Residential 
HAR = High Activity Residential 
HAC = High Activity Commercial 
CBD = Central Business District 
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month. At the end of the year's simulation, an annual total is cal- 
culated . 


6.5 Si mulation of Magothy River with ECHOS Model 

The ECHOS model approximates the natural basins of irregular shape 
as areally equivalent rectilinear areas having areally averaged hydro- 
logic parameters derived from existing data. These segments are tri- 
butary to channel segments with assumed rectangular cross sections and 
uniform geometric and hydrologic characteristics all along its length. 
The numbers of blocks and stream segments required for simulation is 
dictated by the size and geometry of the actual basin. 

Figure 19 shows the Dividing Creek sub-basin. The sub-basin is 
composed of one main stream separating two half basins. Both half 
basins generally flow into Dividing Creek rather than into the Magothy 
River. The sub-basin was approximated by two blocks and a single stream 
segment. Figure 20 shows the computer approximation for the Little 
Magothy River sub-basin. Because a significant fraction off this sub- 
basin drains directly into the main channel of the Magothy River, it 
was approximated by three blocks as shown in the figure. 

Average hydrologic, soils and physiographic properties were as- 
signed to each sub-basin as shown in Table XIII, The computer program 
reads in the parameters for a single sub-basin, the basin structure, 
monthly rainfall profiles, and the loading factors and then performs 
overland flow and channel routing calculations, 

A series of twelve monthly simulations were made for each of the 
sub-basins. The results from each run were converted to pollutant 


FIGURE 19 

COMPUTER REPRESENTATION OF THE 
DIVIDING CREEK WATERSHED 


iCrv5lii' 



ACTUAL 


SIMULATED 




{/' 
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^ " Fir.URE 20 ^ 

COMPUTER REPFESENTATION OF THE 
LITTLE MAQOTHY RIVER WATERSHED 
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XABLE XIII 


MAGOTHY RIVER BASIN - MODEL ECHOS MCCEL DATA INPUT 






OVERLAHO 

FLOW 





INFILTRATICH 


SUBBASIM • 



SLOPE 

{*) 

AVERAGE 

lEHCTH 

(■) 

CHANNEL 

LENCHH 

(«) 

WIDIH 

(-) 

mannings 

"N" 

SOIL 

CMACITY 

(a) 

m/hr 

PEBOOLATIOM 

1 

w 

494423.79 

.074 

338.08 


1507 

.072 

.023 


.001 


E 

483271.37 

.037 

272.07 

* 

1776 

.065 

.064 


.001 

2 


8744909.78 

.017 

540.14 

* 

12487 

.066 

.038 

.065 

.001 

3 

W 

2628253.78 

.038 

486.12 

2743.90 

5407 

.076 

.039 

.054 

.001 


B 

2226765.79 

.017 

606.15 

3674 

.064 

.055 

.069 

.001 

4 

B 

825278.81 

.015 

544.14 

2743.90 

1517 

.083 

.034 

-059 

.001 


B 

1799256.50 

.039 

640.16 

2811 

.098 

.031 

.049 

.001 

S 

B 

3936802.95 

.021 

720.18 

3963.41 

5466 

.071 

.038 

.121 

.001 


fl 

3935687.71 

.008 

558.14 

7051 

.078 

.043 

■ 113 

.001 

6 

B 

11192803.06 

.035 

1104.28 

7621.95 

10136 

.070 

.043 

.140 

.001 


E 

12695419.07 

.021 

849.81 

14939 

.064 

.044 

• IM 

.001 

7 

H 

1286245.35 

.027 

449.75 

3048. 78 

2860 

.078 

.048 

.116 

.001 


E 

1037174.72 

.028 

426.83 

2430 

.093 

.053 

.138 

.001 

B 

H 

407B066.B9 

.025 

943.34 

6097.56 

4328 

.077 

.043 

.121 

.001 


E 

28884 75. 82 

.033 

576.14 

5014 

.068 

.047 

.082 

.001 

9 

W 

2200743.48 

.017 

800.30 

2743.90 

2750 

.052 

.041 

.123 

.001 


E 

817843.86 

.034 

360.09 

2271 

.054 

.052 

.082 

.002 

10 

W 

1148698.88 

.029 

576.15 

2134.15 

1994 

.084 

.033 

.040 

.001 


E 

360 594 . 79 

.055 

X4.08 

1186 

.042 

.038 

.068 

.001 

11 

H 

3713754.63 

.037 

594.15 

3963.41 

4568 

.085 

.044 

.050 

.001 


E 

1237918.31 

.033 

504.13 


2456 

.056 

.050 

.045 

.001 

13 

M 

1858736.05 

.036 

117.30 

4115.85 

15846 

.057 

.029 

.047 

.001 


E 

2788104.08 

.049 

534.13 

5220 

.061 

.046 

.061 

.001 

13 

W 

2174721.18 

.043 

600.15 

274 3.90 

3624 

.073 

.055 

.122 

.001 



2068439.00 

.069 

416.10 

4971 

.081 

.040 

.056 

.001 



1360223.47 

.032 

640.16 

* 

2108 

.039 

.030 

.046 

.001 

14 

H 

3572490.69 

.034 

792.20 

3048.78 

4510 

.073 

.051 

.057 

.001 


E 

3468401.47 

.051 

492.12 

5016 

.042 

.039 

.077 

.001 

IS 

N 

2068439.00 

.026 

368.09 

3963.41 

5619 

.067 

.043 

.076 

.001 



2117152.34 

.013 

688.17 


3077 

.047 

.059 

.109 

.001 


s 

1199095.07 

. 

.006 

496.12 

* 

2417 

.082 

.061 

.090 

.001 




* 

Sub-basln drains directly 

Into KAgothy River 
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loads using the two transfer tables supplied by the RPC. 1-Ionthly totals 
were computed for each sub-basin by multiplying the simulation values 
by the average number of storms for that month and an aggregate yearly 
total was determined. The simulation of the Magothy was performed for 
each of the 17 sub-basins shown delineated in Figure 21. 

6 . 6 Non-Point Source Pollutant Loads for the Magothy for 1977 

Average monthly runoff was simulated for the 1977 calendar year 
for the Magothy Basin in Anne Arundel County, Maryland. 

Pollutant concentration coefficients derived from empirical data 
and LANDSAT derived surface cover were used to convert monthly runoff 
to monthly pollutant load for the 17 sub-basins comprising the Magothy 
Basin. 

Monthly and annual pollutant loads were estimated including 
total suspended solids, BOD, COD, total Kjedehl' nitrogen, total 
nitrogen and total phosphorous using average loading coefficients 
supplied by the RPC. 

Individual storm load estimates compared very closely to the 
limited stream sampling data supplied by the RPC. 

Table XIV summariaes • the . % of- annual pollutant - loads es timated 
by the ECHOS model for the Magothy Basin by sub-basin. The results 
of the simulation indicates that the storm load contribution for the 
pollutants in metric tons /year are given in Table XV. 

Table XVI presents the annual average unit pollutant loads for the 
sub-basins and represents the pollutant loads adjusted to compensate 
for the area variation between the sub-basins. When the unit pollutant 
loads for each pollutant component, is^compare-d to the unit volume of 
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figure 21 

MAGOTHY RIVER SUB-BASiNS 



SUBBASIN SUBBASIN 
#1 NAME 


Otter Pond 
Cornfield Cl 
Grays Ck. 
Blackhole Cl 
Cockey Ck. 
Lake Waterft 
Old Man Ck, 
Cattail Ck. 
Cape Arthur 


AREA 

% 

SUBBASIN 

SUBBASIN 

AREA 

% 

(hectares) 

AREA 

#II 

NAME 

(hectares) 

AREA 

116.2 

1,4 

10 

Cypress Ck. 

151.0 

18 

669.2 

7.8 

11 

Dividing Ck. 

395.0 

4.6 

485.6 

5.6 

12 

Mill Ck. 

464.8 

5.4 

262,6 

3.1 

13 

Forked Ck. 

435.0 

5.1 

757.4 

8.8 

13B 

Bayberry 

122.6 

1.4| 

2370.1 

27.6 

14 

Deep Ck . 

604.1 

6.9 

232. 3 

2. 7 

15 

Little 

Magothy River 

426.3 

4.9 

497.0 

8.1 

15B 

Ty dings 

108.1 

1.3 

302,0 

3.5 
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TABLE XIV 

RPC - TRIANGLE J DATA 
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TABLE XIV (contM) 
WASHINGTON - COG DATA 


% OF TOTAL POLLUTANT LOAD 


AREAS 

TSS 

BOD 

COD 

TKN 

TN 

TP 

AREA 

Otter Pond I 

1.3 

- 

2.7 

3.4 

3.3 

2.3 

1.3 

Cornfield 

5.3 

- 

5.9 

5.9 

5.8 

5.5 

7.8 

Grays Creek 

2.4 

- 

3.2 

4.2 

4.2 

3.1 

5.7 

Blackhole Creek 

,3 

- 

.6 

.75 

.63 

.78 

3.0 

Cockey Creek 

3.2 

1 

4.1 

4.2 

4.1 

3.1 

8.8 

Lake Waterford 

34.2 


25.4 

24.3 

i 

25.6 

1 

27.3 

27.6 

Old Man Creek 

.36 

- 

.8 

1.0 

1.0 

.8 

2.7 

Cattail Creek 

.75 

- 

1.0 

1.1 

1.05 

.8 

8.1 

Cape Arthur 

1.37 

- 

1.4 

1.9 

1.9 

1.6 

3.5 

Cypress Creek 

1.49 

- 

2.5 

3.0 

3.05 

2.3 

1.8 

Dividing Creek 

1.05 

- 

1.8 

1.6 

1.5 

1.6 

4.6 

Mill Creek 

9.8 

- 

10.0 

8.5 

8.2 

8.1 

5.4 

Forked Creek 

6.9 

- 

9.2 

8.5 

7.9 

7.0 

5.1 

Bayberry 

8.2 

- 


7.2 

7.5 

7.0 

7.0 

Deep Creek 

16.9 

- 


18.6 

18.7 

16.0 

5.0 

Little Magothy River 

5.7 

- 


3.6 

3.6 

4.7 

1.4 

Tydings 

.9 

- 

1.8 

2.2 

2.1 

1.6 

1.2 

MT/yr 


|H| 

498.8 

8.02 

9.50 

1,28 


Kg/ha. 


n 

58.0 

0.93 

1.11 

0.15 





























TABLE XVI 


AVERAGE ANNUAL POLLUTION LOADING (Kg/Ha) MAGOTHY RIVER BASIN 



SUB-BASIN 

TSS 

BOD 

COD 

TKN 

TN 

TP 

volume 

DISaiMGE 

UNIT 
VOLUME 
(m /hec) 

1. 

Otter Pond 

164.82 

9.88 

91.51 

1.95 

2.66 

0.43 

19A29.06 

167.20 

2. 

Cornfield 

51.39 

2.28 

30. 32 

0.56 

.82 

0.12 

23600.95 

35.27 

3. 

Grays Creek 

50.97 

1,68 

26.34 

0.60 

.84 

0.12 

23236.03 

43.73 

4. 

Black Hole Creek 

14.28 

1.04 

8.76 

0.17 

.22 

0.04 

4267.46 

16.27 

5. 

Cockey Creek 

32.40 

3.24 

11.20 

0.72 

1.02 

0.08 

20444.90 

26.99 

6. 

Lake Waterford 

70.54 

2.58 

35.83 

0,63 

1.02 

0.15 

105387.43 

44.47 

7. 

Old Man Creek 

24.88 

1.58 

14.14 

0.29 

.40 

0.07 

6065.37 

26.11 

8. 

Cattail Creek 

8.99 

.92 

5.20 

0.10 

.28 

0.02 

5088.14 

7.30 

9. 

Cape Arthur 

38.54 

1,74 

18.82 

0.43 

.60 

0.08 

8563.79 

28.36 

10. 

Cypress Creek 

121.31 

6.96 

65.89 

1.37 

1.90 

0.30 

17322.87 

114.72 

11. 

Dividing Creek 

22.38 

1.64 

15.24 

0.24 

.68 

0.06 

8444.76 

21.39 

12. 

Mill Creek 

108.81 

4.16 

68.52 

1.11 

1.68 

0.24 

31231.09 

67.19 

13. 

Forked Creek 

76.09 

5.64 

71.09 

1.18 

1.72 

0.28 

41534.26 

95.48 

13B 

Bayberry 

S60 .'30 

12.15 

195.38 

3.89 

5.82 

0.73 

2375.25 

193.78 

14. 

Deep Creek 

186.21 

8.04 

104.55 

2.02 

2.94 

0.42 

80810.68 

133.77 

15. 

Little Magothy 
River 

52.21 

1.56 

32.33 

0.47 

.80 

0.11 

10953.92 

25.70 

15B 

Tydings 

116.83 

7.12 

65.03' 

1.38 

1.88 

0.30 

13371.38 

123.70 
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runoff, it is evident that the primary generator of high unit pollu- 
tant load response is the hydrologic response of the basin, i.e., the 
quantity of runoff per unit area. In each case, except Mill Creek, 
there is a direct correspondence between high unit pollutant load and 
high unit runoff. 

The relatively high pollutant loads from Mill Creek can possibly 

be explained by the high fraction of pasture and commercial land use 

in Mill Creek which contribute high total suspended solids. Contrast 

the pollutant loads from Mill Creek with a unit volume of runoff of 
3 

67.2 m /ha. with that of the Tydings sub-basin which has equivalent 

3 

pollutant loads but has almost double the (123.7 m /ha) unit volume 
of runoff. 

Figure 22 displays the simulated monthly runoff for selected sub- 
basins of the Magothy. The relative magnitude of the curves are a 
direct result of the area variation since the rain was considered to 
be uniform over the total area. These curves all have a similar shape 
which is dictated by the assumption of uniform rainfall rate throughout 
the Magothy Basin, Minor differences, particularly for low values of 
runoff are a result of the variation of soils and infiltration between 
the sub-basins. For example, Blackhole Creek, which has approximately 
3% of the total area,, has significantly less runoff for the period 
January through June. Blackhole Creek is comprised of soils with 
higher percolation rates. Cockey Creek* which is approximately 9% of 
the total area, e^diibits relatively high runoff for the whole year. 



MOmXY DISCHARGE, m5(xlO^ 
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FIGURE 22 

COMPUTED MONTHLY RUNOFF 

FDR MA30THY RIVER SUB- BASINS WITH FLOW-DOMINATED 

POLLUTANT LOADING 



TIKE .MDNTHS 


LEGEND ; 

gray's creek 

BLACKHOLE CREEK 

COCKEY CREEK 

OLDMAN CREEK 

CATAIL CREEK 

— ■» — BAYBERRY BASIN 
TYDINGS BASIN 
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Since the pollutant load conversions are sensitive to the volume 
of runoff, the monthly runoff curves are a good indicator of the timing 
of the major pollutant contributions. 

Figure 23 shows the relative monthly contribution of pollutants 
for the Dividing Creek sub-basin and is characteristic of all of the 
sub-basins. The large contribution for the period centered upon 
August is a direct consequence of the rain profiles characteristic of 
the post summer period. Table XVII presents the summary of the average 
monthly rainfall event profiles. The table shows that the March-May 
period is characterized by low intensity rains of relatively long dura- 
tion and which have relative low unit volumes. This combination tends 
to produce low runoff and hence low pollutant loading contributions. 

The July-September period, in contrast, is characterized by relatively 
short, high intensity rain events which produce higher runoff and thereby 
contribute high pollutant loads. The shift of the peak period toward 
the winter months is probably due to the seasonal effect of evapotrans- 
piration which peaks in July and rapidly decays toward December as 
shown- in Figure 24. 

Figures 25a through f show the annual unit pollutant load (kg/ 
hectare) for the Magothy Basins for the pollutants considered. 

These figures are compared with similar ones for land use and runoff in 
Figures 26 to 29. It is noted that the values for pollutant loads are 
rationally explained for the most part by the values of unit runoff. 

High values of unit runoff yield high across-the-board pollutant 
loadings j and the reverse is true for low unit runoff values . In cer- 
tain cases, however, pollutant loading is "driven” by land use rather 



FIGURE 23 


TIME PHASING OF POLLUTANT LOAD 



MONTHS 



Cumulative fraction of 
total pollutant load 


Monthly fraction of 
total pollutant load 





January 

DURATION 

9.12 

February 

8.96 

March 

7.55 

April 

6.40 

May 

4,93 

June 

3.52 

July 

3.32 

August 

4.18 

September 

5.39 

October 

7.62 

November 

8.38 

December 

10.91 
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ABLE XVII 


ihe average rainfall event 


UNIT 

INTENSITY 

VOLUME 

mm/hr. 

7.9 

.9 

9.8 

1.1 

9.9 

1.3 

7.8 

1.3 

8.8 

1.8 

9.5 

2.7 

10.8 

3.2 

12.3 

2.9 

10.7 

2.0 

11.1 

1.5 

10.9 

1.3 


12.3 




EVAPORTRANSHRATION, INCHES 
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UNIT POLLUTANT LOADING - MAGOTHY RIVER SUBBASINS 
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UNIT POLLUTANT LOADING - MAGOTHY RIVER SUBBASINS 
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unit POLLUTANT LOADING - MAGOTHY RIVER SUBBASINS 
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FIGURE 27 


PERCENT URBAN /RESIDENTIAL AREA - MAGOTHY RIVER SUB-BASIN 


URBAN-RESIDENTIAL 
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FIGURE 28 

PERCENT AGRICULTURE - PASTURE AREA - MAGOTHY RIVER SUBBASINS 
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FIGURE 29 
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than flow. As shown in Table XVIII, seven of the Magothy River sub-basins 
produced pollutant loads which could not fully be explained by unit dis- 
charge. Mill Creek, for example, did not behave like the flow dominant 
Blackhole sub-basin. The unexpectedly high values of suspended solids, 
COD, TKN, and phosphorous result from the high proportion of pasture 
present rather than unit discharge. This exemplifies the importance of 
adequate land use data to non-point source pollutant modelling. 



TABLE 


MAGOTHY RIVER SUB-BAS INS - DRI 


SUBBASIN 

TSS 

BOD 

COD 

1. 

Otter Pond 

H 1 

H 

H 

2. 

Cornfield 

H-M 

M 

M 

3. 

Gray's Creek 

H-M 

M 

M 

4. 

Black Hole Creek 

L 

L 

L 

5. 

Cockey Creek 

M 

M 

L 

6. 

Lake Waterford 

H 

M 

M 

7. 

Old Man Creek 

L-M 

L 

L 

8. 

Cattail Creek 

L 

L 

L 

9. 

Cape Arthur 

M 

L 

M 

10. 

Cypress Creek 

H 

H 

H 

11. 

Dividing Creek 

L-M 

L-M 

L-M 

12. 

Mill Creek 

H 

M 

H 

13. 

Forked Creek 

H 

M 

H 

13B. 

Bayberry 

H 

H 

H 

14. 

Deep Creek 

H 

H 

H 

15. 

L. Magothy River 

H 

L 

M 

15B. 

Ty dings 

H 

H 

H 







III 


VERS OF POLLUTANT LOADING 




UNIT 

POLLUTANT LOAD 

TN 


DISCHARGE 

DRIVER 

H 

H 

H 

Runoff 

L-M 

M 

M 

Runoff 

L-M 

M 

M 

Runoff 

L 

■ 

L 

Pvunof f 

M 

B 

L 

Medium Values due to 7% 
Pasture 

M 

M 

M 

High TSS due to 6% agri- 
culture 

L 

L 

L 

Runoff 

L 

L 

L 

Runoff 

L 

L 

L 

Medium values due to 
12% strip commercial 

H-M 

H 

H 

Runoff 

L 

L 

L 

Low-Medium values due to 
1% strip commercial 

M 

H 

M 

High values due to 20% 
pasture 

M 

H 

M 

High values due to 11% 
pasture, 3% strip com- 
mercial 

H 

H 

H 

Runoff 

H 

H 

H 

Runoff 

L 

M 

L 

Medium, High values due to 
6% agriculture; 25% pasture 

M-H 

H 

H 

Runoff 


I 

c» 

OJ 

I 
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VII. THE USE OF THE LANDSAT DERIVED POLLUTANT LOAD AT EPC 

The resolution of region-wide land cover measurement and concom- 
mitant non-point source pollutant load assessment has had a beneficial 
impact upon "208" planning work. Prior to demonstration of optical 
classification techniques using LANDSAT data, a major concern over pro- 
ject implementation had existed. The RPC now feels it has at its disposal 
techniques which greatly reduce the analytic effort while maintaining 
satisfactory precision. 

Using a hybrid procedure of LANDSAT and aerial remote sensing, 
land use analysis can be accomplished to the required level of detail 
(Level 1) in very short time for small areas. Basins of size compar- 
able to the Magothy could be analyzed by one person in approximately 
three days after the imagery has been received and prepared. An op- 
portunity for resource savings exists. Using computer classification 
of LANDSAT CCT, accurate land use data for the total region and larger 
areas within the region have been generated. 

Response to these procedures by the RPC personnel has been favor- 
able. The work performed on the Magothy has been included in the RPC’s 
draft "208” plan as indicated in the following section. 

7 . 1 Excerpt from the -RPC "208" Plan 

Magothy River Pollutant Load Analysis (Note figures not included) 

The analysis of water quality samples during wet weather condi- 
tions is inportant for an accurate understanding of non-point sources 
of pollution. The cost of wet weather quality sampling, however, is 


extremely high. 




An alternative to the actual collection of wet weather water 


quality data is the development and calibration of a computer model 
which simulates natural stream conditions. This approach is being 
demonstrated through a special analysis of the Magothy River Watershed. 

The computer model developed for the Magothy River Watershed uses 
land cover data from LANDSAT satellite imagery and aerial photography, 
pollutant load data, and rainfall and runoff records. The model cal- 
culates wet weather pollutant loadings, using these data. 

The steps followed in developing and applying the computer model 
for the Magothy River are listed below: 

1 . Land Use Analysis from LANDSAT and Aerial Imagery : 

The remote sensing analysis of the Magothy River watershed, was 
completed during the January reporting period. 

The procedure first required the preparation of satellite and 
aerial imagery into color slide form, as described in the 
December report. Subsequently, a preliminary examination of 
the aerial scene containing the Magothy River was made. It 
was determined that seven general land use classes were present: 
residential, commercial/industrial, extractive, agriculture, 
water, forests, and pasture/grass areas. 

The entire Magothy River Watershed was mapped and classified 
from the aerial photo. Because of the age of the aerials 
(September, 1970), Imagery from 1975 was acquired from the 
Baltimore RPC, and the ground truth map updated. Changes had 
to be made in about ten areas reflecting new residential or com- 
mercial development. 

Test sites were then selected for each class. This area con- 
stituted about two to three percent of the total watershed. 

A map was prepared locating the sites and the watershed boun- 
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This map was used as the basis for LANDSAT classifications of 
land use. Each satellite scene was projected over the map of 
test sites and the colors of each class were recorded. The 
images of May, August, and October, 1975, February, 1974, and 
April and October, 1973 were used. This information was then 
used to determine which scenes were preferable for identifying 
different land use classes. For example, forest appeared dark 
in the February, 1974, scene while potential confusers were 
lighter. This image was therefore used to segregate forest. 
Similar logic was applied to separation of the other land use 
classes , 

High color contrast in the October 6, 1973 scene allowed the 
relatively easy separation of the extractive and commercial/ 
industrial classes. The negative image of October 6 showed 
water as bright white; it was therefore readily separated. 

Three classes - residential, pas ture/ grass , and agricultural - 
then remained to be segregated. No single image shows agri- 
cultural areas in marked contrast with the other classes. 
Therefore, the temporal characteristics of cropped fields 
had to be relied upon. By examination of all the available 
scenes, it was determined that agricultural areas followed a 
unique progression of colors: dark pink in April, light blue 

to white in October, and orange/light blue in February. All 
areas exhibiting this behavior were then classified as agri- 
cultural. No scene indicates residential areas strongly 
separable from pasture/grass. This is due to the low density 
(2-4 units per acre) of most residential development present. 
Some of the more densely developed residential areas could be 
consistently identified, but most behaved like pasture/grass 
areas all through the year. The total residential fraction 
could not be segregated to the accuracies required by the 
Regional Planning Council. It was decided, therefore, to ex- 
tract the residential class from the aerial photography. While 
the use of computer compatible tapes or additional images might 
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permit the total classification from LANDS AT, the use of a 
hybrid satellite - aerial remote sensing system to derive 
pollution loads represents a tool useful to the Regional 
Planning Council. 

On the basis of this technique, the land use map shown in 
Figure 3-B-23 was prepared. Overall classification error 
was 2.3%. 

2 . Pollution Loading Factors 

Two sets of pollution loading figures received from the RPC 
and the Washington COG were used in the modeling. These sets 
of numbers were reasonably consistent; they were used as 
multipliers to convert flow to kilograms of pollutants. 

3. Rainfall/Runoff Records 

Runoff records were received from the RPC for three ‘separate 
sites: Cape St. Clair, Bay Hills, and Sevema Park. Suf- 

ficient data points to plot outflow hydrographs, which 
matched available rainfall data, existed only for four events 
at Cape St. Clair. 

These events were for February 24, March 18, April 2, and 
July 25, 1977. 

In addition to on-site rainfall data, records were obtained 
for the Baltimore gage, giving hourly rainfall values for 
1948-1974. From these data, average monthly rainfall pro- 
fits were developed. 

4. Model Calibration 

The Ecosystems Hydrologic Simulator (ECHOS) was calibrated 
using the data from Cape St. Clair. The physical parameters 
of area, slope, flow length, and stream characteristics were 
taken from USGS 1:24000 topographic maps. Land cover (sur- 
face friction) came directly from the map prepared from the 
satellite and aerial analysis. 
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The four rainfall events were run and matched to actual 
hydrographs. Runs were made until a satisfactory fit was 
achieved. An example is shown in Table 3-B-15. These runs 
were used primarily to set the values of subsurface soil 
moisture depletion. 

The four events, combined with the Baltimore rainfall data, 
were also used to set appropriate antecedent soil moisture. 
From the discharge record, the percentage of rainfall which 
became runoff was computed. The fractions derived are shown 
below : 


February 

1.0% 

March 

2.0% 

April 

1.3% 

July 

2.8% 


From these four points, curves were extrapolated based upon 
comparison of rainfall columns and peak rates as derived from 
the 25 year Baltimore records. 

It was now known, for example, that in August about 2.8 percent 
of rainfall would become direct runoff on the average. A second 
set of calibrations was then performed to set antecedent soil 
moisture (API) to a value that would produce the correct per- 
centage of runoff. This was done by trial and error, setting 
an estimated value of API and recording the resultant runoff. 

The optimal set of antecedent soil moisture percent are given 
below: 


Month 

Optimal % API 

January 

99% 

February 

99 

March 

100 

April 

100 

May 

96 

June 

94 

July 

93 

August 

93 
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September 

95 

October 

96 

November 

97 

December 

98 


The values derived are very well balanced when considered from 
a hydrologic point of view. High values of soil moisture were 
obtained for months with lower solar radiation and ET and vice- 
versa. 

These values were input to the ECHOS data files. The simula- 
tion of pollutant loads could now be run. 

5 . Model Runs 

The Magothy Basin was broken down into major sub-basins. The 
initial run was made for sub-basin 11, Dividing Creek. ECHOS 
reads the rainfall profile month by month, calculates overland 
flow, routes the runoff through the stream channel and computes 
the hydrograph at the outlet. Streamflow is cumulated and con- 
verted to pollutant loads using the COG and RPC factors des- 
cribed earlier. Monthly loads and an average annual summary 
are produced* 

Total loads are given by multiplying the load from the average 
storm by the average number of storms per month as calculated 
from the Baltimore data. An example of the output from the 
computer analysis runoff summary and pollutant loadings for 
the Dividing Basin during July is shown on Table 3-B-16. 

7.2 The RPC^s Direct Response 


The favorable reaction to the Magothy River analysis stems directly 
from the fact that the technique was demonstrated for an area where the 
results would be immediately useful. 

The written comments given in Appendix B were received from the RPC on 
the use of LANDSAT data to derive non-point source pollutant loads. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


The application of LANDSAT derived data to non-point source pol- 
lutant load estimation of the Magothy Basin was successfully accom- 
plished for the Baltimore Regional Planning Council. The pollutant 
loads and the LANDSAT- oriented techniques which were used to derive 
them have been incorporated into the current RPC "208” plan and have 
been selected for future broader applications at RPC. The cost-ef- 
fectiveness and the applicability of the techniques for present and 
future use in regional NPS investigations have been judged to be ex- 
cellent by the RPC personnel. 

A. Significant Findings as Regards LANDSAT Derived Land Cover: 

1. The relatively high cost of aerial photography and the 
limited frequency of total regional coverage gave rise 
to the need for LANDSAT-oriented land cover classifica- 
tion techniques. The need for accurate sub-basin land 
use data requires an improvement over conventional land 
use estimates with large cell size. LANDSAT computer 
classifications are adequate for large area estimates, but 
exhibit relatively high errors for small areas. The de- 
velopment of a hybrid, LANDS AT-opti cal technique with 
minimum aerial coverage input for residential areas 
proved a successful, cost-effective solution to the pro- 
blems in the Baltimore RPC area. 

2. Average errors of less than 5% were obtained for the 
four small areas investigated in these studies. The 
area evaluated ranged from 108 hectares to 8,600 hec- 
tares . 

3. Accuracies obtained by this hybrid analysis technique 
were commensurate with the error budget adopted for non- 
point source pollutant load estimation by the RPC. 

B. Significant Findings as Regards Pollutant Loads for the 

Magothy Protypical Area Include: 

1, Annual pollutant loads from storms resulting from the 
average rainfall were computed to be: 615.3 metric 

tons per year of total suspended solids, 26.5 metric 
tons per year of BOD, 344.7 metric tons per year of 
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COD, 6.3 metric tons per year of TKN, 9.5 metric tons 
per year of total nitrogen and 1.4 metric tons per year 
of total phosphorous. 

2. The major period of storm load contribution was the 
period of July through September. The characteristic 
occurrence of high intensity storms in this period ac- 
counted for 70% of the total storm load. 

3. High unit pollutant loads were computed for the following 

sub-basins of the Magothy: Otter Pond, Cypress Creek, 

Mill Creek, Bayberry, Deep Creek and Tydings. These sub- 
basins total 31% of the total area of the Magothy. 

C. Significant Findings as Regards the Application of Hydrologic 

Models Based Upon Satellite Data for Non-Point Source Load 

Estimation 

1. A simple, useful, and cost-effective modeling procedure 
has been developed to estimate annual and event-oriented 
pollutant loads from storms, using specific or generalized 
pollutant load coefficients. 

2. The procedure is based upon land cover data derived from 
satellite images. 

3. The accuracy levels of the procedure is commensurate with 
the accuracy level of the input pollutant data. 

4. The procedure provides environmental and land use 
planners with practical, easily adopted techniques 
from which to estimate the sensitivity of specific 
watersheds to current and future land use distributions 
and patterns, 

5. The model developed has been made extremely flexible to 
accommodate the maximum user-derived inputs, 

6. Techniques for estimating the input parameters for areas 
with limited data have also been developed and demon- 
strated. 

The experience gained during this technology transfer project 
supports the expectation which surround the availability of frequent 
accurate land cover data from satellites such as LANDSAT. The user 
population is large and widespread. Certain problems like those en- 
gendered by the EPA "208" program are immediately amenable to solutions 
using LANDSAT data. 

k. 
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Technology transfer requires that not only the user be made 
familiar with the LANDS AT products but that he also be aided in de- 
veloping the models, techniques and algorithms necessary for its 
successful employment. Most users are not engaged in pure research 
pursuits and hence require "proof" before committing these resources 
to the use of a new and not as yet proven sources of data. Once 
these problems are overcome, however, the users are quick to grasp and 
improve the use of the data. 

1. It is recommended that the procedures developed herein be made 
available to a broad spectrum of "20R*'-type users. 

2. It is further recommended that the utility of pollutant load 
estimation techniques be tested for the so-called non-designated 
areas throughout the U.S. where data is much less available than 
in designated regional areas such as the Baltimore EPC. 

3. The correlative parameter in this pollutant load model was 
chosen to be total runoff from an event due to the unavail- 
ability of flow to load correlations. It is recommended 
that further simulations be carried on as these data become 
available and that a concomitant expansion of the rainfall 
recurrence case be considered to include expected non-lin- 
earities in the flow-to-load relationships. 

A. Following the suggestions given by the Baltimore RPC, it is 
recommended that a procedural handbook be prepared illus- 
trating the application of LANDSAT data to non-point pol- 
lution estimation for subsequent transmittal to the large 
user audience throughout the U.S. 



APPENDIX A: ECHOS MODEL DESCRIPTION 
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Description of the' Model 

The Ecosystems Hydrologic Simulator (ECHOS) requires 
that the following data about each sub-basin be known or 
estimated: 

1. Average slope, m/m 

2 

2. Area, m 

3. Average Overland Flow Length, m 

4. Width (equal to area v flow length), m 

5. Average Mannings "n" 

6- Soil Moisture Capacity, m 

7. Antecedent Moisture Content, m 

8. Final (saturated) Infiltration Rate, m/sec. 

9. Percolation Rate, m/m 

10. Detention Storage, m 

Additionally, these data are required for each stream 
s egemen t; 

1. Length, m 

2. Average slope, m 

3. Channel width, m 

4. Channel Manning's "n" 

As adapted for the Baltimore RPC, (Version IP) data 
are also required giving pollutant loading factors for two 
sets of land use classification, one supplied by RPC and 
One received from the Washington Council of Government (COG) 



RPC Classes 


Rural (^4% impervious) 

Low Activity Commercial 12% impervious) 

Low Activity Residential 16% impervious) 

High Activity Commercial 35% impervious) 

High Activity Residential 32% impervious) 

Central Business District 80% impervious) 

COG Classes 
Shopping Centers 
Strip Commercial 
High Rise Residential 
Town House/ Garden Apartments 
Med. Density Single Family Residential 
Construction 
Forest 

Conventional Tillage Agriculture 
Minimum Tillage Agriculture 
Pasture 

The program is written in language. It is presently 


configured for a PDP 11-34, 32K computer but has also run 

on IBM and Honeywell machines of comparable size. 



The program reads th'® data it requires out of separate 
files prepared in advance of a run on disk or tape. Four 
files are involved; for the Magothy River they have been 
called DATA.DAT, RAIN.DAT, LINK.DAT, and CON. DAT. Examples 
of each are given in Table AI a,b,c, and d. 

DATA.DAT contains values for all the physiographic 
parameters of each sub-basin and stream. It also includes the 
initial conditions of the simulation and time-keeping informa- 
tion. RAIN.DAT contains a set of intensities and corresponding 
durations for each month. LINK.DAT gives the program the 
connection between sub-basins and channels, i.e., it identifies 
which basins flow into which streams. CON. DAT contains the 
pollutant loading factors and land use fractions for the RPC 
and COG data. 

A simplified flow chart for ECHOS is shown in Figure lA. 
Briefly, the model reads the first rainfall profile and the 
parameters for the first sub-basin. It computes overland 
flow and stores it as elements of a f low- ve r sus- t ime matrix. 

The process is repeated for each sub-basin. Subsequently, the 
model rates the output of the appropriate sub-basin or sub-basins 
through the channel or channels as specified by LINK.DAT. 

Flow is routed until the basin outlet has been reached. A 
table is then printed out giving time since the start of 
rainfall, rain and discharge rate, and cumulative rain and 
discharge. Finally, total discharges for the month are computed 
and converted to pollutant loads using the loading factors 


TABLE A1 


DATA FILES 


REENTER 

*ERDX1 :C0N.DAT*EUDX1 :CON.DAT»RiBi/Li» 

04191006.300000000.0025:13903.300000000.000000000,000000000.0 

0167.000007.000038.200000.710000.41 

0239.000007.600061 . 100001 .030000.47 

0293 . 000007 . 700076 ,100002 . 320000 . 49 

0751.000012.200140.600002.340000.68 

0633.000013.200172.800002.220000.60 

0179,000013. 300178.200002 .970000.46 

00000000,000000000.000000000.000000000.000000000.002553903 . 3 

00000000.003443794.600000000.000000000.000747211 .9 

0074.000054.660001.810002.430000.20 

0098.000111.530001 .330001 .760000,27 

0034.000096.280000,870001 .200000.22 

0115.000073.460001 . 120001 .490000.29 

0042.000095.880001 .940002.420000.38 

0142.000115.150002,760003.640000,35 

0763.000089.350001.610001.760000.26 

0000.000072.400001 .280001.300000.12 

1681,000853.380008.900013.420003.18 

0133.000086.240001.930003.520001.14 

0544.000549.000004.970003.160001 .06 


.R EPIT 

♦ ERDXl :DATA.DAT»EUDX1 iDATA.DA-CtRSBf/L** 

01015.0600000000.0. 00000 

00540.140.0170.066012478. . 002500000.0380000000.0000002420.000018056 
CORNFUO 06744909.8 

01 

01000.000. 0040.0382.00 
*EX«$ 


R EDIT 

tERDXi;LlNK.DAT»EUDXl!LINK.DAT4R$B$/L** 

EAST 

OOOIDUMMY 

01 

4EX«$ 


A 

00.9908. 14JANUARY 18032400 

.00000003601800. 

.00000010701800. 

.00000017901800. 

.00000025001800. 

.00000032201800. 

,00000039301800. 

.00000047201800. 

.00000045101800. 

.00000041501800, 

.00000037201800. 

.00000032901800. 

.00000028601800. 

.00000023601800. 

.00000019301800. 

,00000015001800. 

.00000010701800. 

.00000006401800. 

.00000002101800. 

00. 9907.63FEHRUARY 16028800 

.00000004201800. 

.00000012501800. 

.00000022201800. 

.00000030601800. 

.00000033401800, 

.00000047201800. 

.00000055601800, 

.00000066701800. 

.00000061101800. 

.00000052001800. 

,00000044401800. 

,00000036101800. 

.00000027001800, 

.00000019401800, 

.00000011101800. 

.00000004201600. 

01 .000B.78MARCH 16020800 

.00000005201800. 
.00000015701800, 
.00000026201800. 
,00000036701800. 
.00000047101800. 
,00000057601800. 
.00000063901800. 
.00000062101800. 
.00000054601800. 
.00000047101800. 
.00000040401800. 
.00000033701800. 
.00000026201800. 
.00000018701800. 
.00000011201800, 
.00000003701800. 
01.0010.24AFRJL 13023400 

.00000005101800. 
.00000016601800. 
.00000027501800. 
.00000039001800. 
.00000049801800. 
.00000060701800. 

, 00000072201UOO. 
.00000052401800. 
.00000042801000. 
.00000033201800, 
.000000236011100. 
.00000014701800, 
.OoOOOOO'JlOlUOO. 


B 

00.9610.52MAY 10018000 

.00000009401800. 
,00000028101000. 
.00000046901800. 
.00000065701800. 
.00000083101000. 
.00000100001800. 
.00000071001000. 
.00000049601000. 
.00000029501800. 
.00000009401800. 

00. 9409. 65 JUNE 14012600 

.00000011000900. 

.00000032900900. 

.00000054900900. 

,00000076800900. 

,00000098800900. 

.00000120700900. 

.00000132000900. 

,00000132000900. 

.00000114400900, 

,00000095600900. 

.00000073700900. 

.00000051700900. 

.00000032900900. 

.00000011000900. 

00.9309, IIJULY 14012600 

.00000012300900. 

,00000035400900. 

.00000058500900. 

.00000081600900. 

.00000104000900. 

.00000127900900. 

.00000151000900 . 

,00000169400<’00. 

.00000132500900. 

,000001 io':oo 9 oo. 

.00000084300900. 

.00000059800900. 

.00000036700900. 

.00000012300900. 

00.9309.46AUGUST 16014400 

.00000010800900. 

.00000032300900. 

.00000053900900, 

,00000077000900. 

,00000100100900. 

.00000121700900. 

.00000143200900. 

.00000159500900. 

.00000159500900. 

.00000146300900. 

.00000123200900. 

.00000100100900. 

.00000078500900. 

.00000055J0090'). 

.00000033900900. 

.00000010800900. 


C 

00. 9507. 29SEFTEm<R 1101 9000 
.00000010401800. 
.00000031201000. 
,00000050501800. 
,00000071301000. 
.00000092101800. 
.00000099001800. 
.00000083201800. 
.00000065301800. 
,00000046001800. 
.00000028201800. 
.00000008901800. 
00.9606.600CTOBER 15027000 
.00000006201800. 
.00000018601800. 
,00000030901000. 
,00000043301800. 
.00000055701000. 
,00000068101800. 
.00000080001800. 
.00000073001800. 
.00000059401800, 
.00000050701800. 
.00000040801800. 
.00000032201800. 
.00000023501800. 
.00000013601800. 
.00000005001800. 

00 . 9706 . 78N0UEMBER 17030600 

.00000005501800. 

,00000014301800. 

.00000023101000. 

.00000033001000. 

.00000042901800. 

.00000051701000. 

.00000061601800. 

.00000085301800. 

.00000059401000. 

,00000051701800. 

,00000045101800. 

.00000038501000. 

.00000031901000. 

.00000024201800. 

.00000017601000. 

.00000009901000. 

.00000003301800. 

00 .9007 . 46IIECEMBER22039600 
.00000003801800. 
.00000011401800. 
.00000018201800. 
.00000025201800. 
.00000032501800. 
.00000040101000. 
.00000047701800. 
.00000104501300. 
.00000054501800. 
.00000050701800. 
.00000046101800. 
,00000042401800. 
.00000038601800. 
.00000034001800. 
.00000030301830, 
.00000026501800. 
,00000021901t‘0(’. 
.00000018201800. 
.00000014401800. 
.00000037801000. 
.00000006101800. 
.00000002301800. 
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FIGURE lA 

ECHOS MODEL; SIMPLIFIED FLOW CHART 











of CON. DAT. Tables showing the monthly loads for organic 
pollutant are then output. (A detailed description of the d 
derivation of ECHOS is given in the final report of NAS8-30539, 
24 February/ 1977) • 

Table All is a sample output of the model for the 
Dividing Creek sub-basin of the Magothy River. The output 
derived from the file data is shown in Table Al a through d. 


Table AIII gives a listing of ECHOS, Version IP. 



TABLE All: SAMPLE OUTPUT; ECHOS MODEL 


ECOSYSTEMS HYDROLOGIC SIMULATOR (ECHOS) 
VERSION IF' 


X subbasin; east * 

X LENGTH* 504.13 METERS SLOPE* 0.033 N= 0.05A UIDTH= 245A. METERS DETENTI0N=0. 00250 METERS * 
X SOIL MOISTURE CAPACITY = 0.050000001 METERS ANTECEDENT SOIL MOISTURE- 0.049500000 METERS X 
X PERCOLATION RATE* 0.000000220 M/SEC FINAL INFILTRATION RATE* 0.000012500 M/SEC X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


X subbasin: uest 

X LENGTH* 594,15 METERS SLOPE* 0.037 N= 0.005 
X SOIL MOISTURE CAPACITY - 0.044000000 METERS 
X PERCOLATION RATE- 0.000000238 M/SEC 


WIDTH* 4560. METERS DETENTION=0 . 00250 METERS 
ANTECEDENT SOIL MOISTURE- 0.043560002 METERS 
FINAL INFILTRATION RATE* 0.000013900 M/SEC 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X STREAM segment: DIVID CK * 
X LENGTH- 3963.41METERS SLOPE* 0.010 N- 0.038 WIDTH- 2.00METERS X 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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++++++++++++++++++++++++++++++++++++++f+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
+ + 
4 4 


4 



RUNOFF 

SUMMARY FOR DIUIDECK BASIN 


+ 

4 



AUG 

. EUENT for: JANUARY 


+ 

4 

4 



AUG. 

NO. OF storms; 

8.14 


+ 

+ 

4 

4 

TIME 

DISCHARGE 


RAIN RATE 

CUMULATIUE 

CUMULATIUE 

+ 

4 

4 

4 

(SECS) 

(M3/S) 


(M/S) 

DISCHARGE(M3) 

RAINFALL(M3) 

+ 

+ 

4 

4 

900. 

0.00000000 


0.0000000360 

0.0000 

128.0342 

+ 

+ 

4 

IQOO. 

0.00000000 


0.0000000360 

0,0000 

256.0684 

+ 

4 

2700. 

0.00000000 


0,0000001070 

0,0000 

636,6145 

+ 

4 

3600. 

0.00000000 


0.0000001070 

0.0000 

1017.1606 

4 

4 

4500. 

0.00000000 


0.0000001790 

0.0000 

1653.7751 

+ 

4 

5400. 

0.00000000 


0.0000001790 

0.0000 

2290.3896 

4- 

4 

6300. 

0.00000000 


0,0000002500 

0.0000 

3179.5161 

+ 

4 

7200. 

0.00000000 


0.0000002500 

0,0000 

4068.6426 

4- 

4 

BIOO. 

0.00000000 


0.0000003220 

0,0000 

5213.8374 

4- 

4 

9000. 

0.00000000 


0.0000003220 

0.0000 

6359.0322 

4- 

4 

9900. 

0.00000000 


0.0000003930 

0,0000 

7756.7388 

4- 

4 

10800. 

0.00000000 


0.0000003930 

0,0000 

9154.4453 

4- 

4 

11700. 

0.00000000 


0.0000004720 

0.0000 

10833.1162 

4- 

4 

12600. 

0.00000000 


0,0000004720 

0.0000 

12511.7871 

4- 

4 

13500. 

0.00000000 


0.0000004510 

0.0000 

14115.7715 

4- 

4 

14400. 

0.00000000 


0,0000004510 

0.0000 

15719.7549 

4- 

4 

15300. 

0.00000000 


0. 0000004150 

0,0000 

17195.7051 

4- 

4 

16200, 

0.00000000 


0.0000004150 

0.0000 

18671.6543 

4- 

4 

17100. 

0.00000000 


0.0000003720 

0.0000 

19994.6738 

4- 

4 

18000. 

0.00000000 


0.0000003720 

0,0000 

21317.6953 

4- 

4 

18900. 

0.00000000 


0.0000003290 

0.0000 

22487.7852 

4- 

4 

19800. 

0,00000000 


0.0000003290 

0.0000 

23657.8770 

+ 

4 

20700. 

0.00000000 


0.0000002860 

0.0000 

24675.0371 

4- 

4 

21600, 

0.00000000 


0.0000002860 

0,0000 

25692.1973 

4 

4 

22500. 

0.00000000 


0.0000002360 

0.0000 

26531.5332 

4 

4 

23400. 

0.00000000 


0,0000002360 

0.0000 

27370.8672 

4 

4 

24300. 

0.00000000 


0.0000001930 

0.0000 

28057.2734 

4 

4 

25200. 

0.00000000 


0.0000001930 

0,0000 

28743.6777 

4 

4 

26100. 

0.00000000 


0.0000001500 

0.0000 

29277.1543 

4 

4 

27000. 

0.00000000 


0.0000001500 

0,0000 

29810.6289 

4 

4 

27900. 

0.00000000 


0.0000001070 

0.0000 

30191.1758 

4 

4 

28800. 

0.00000000 


0,0000001070 

0.0000 

30571.7207 

4 

4 

29700. 

0.00000000 


0.0000000640 

0.0000 

30799.3379 

4 

4 

30600. 

0.00000000 


0.0000000640 

0.0000 

31026.9531 

4 

4 

31500. 

0.00000000 


0.0000000210 

0.0000 

31101.6406 

4 

4 

32400. 

0.00000000 


0.0000000210 

0.0000 

31176.3262 

4 


+ 4 
4 4 
4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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++++++++++++++++++++++♦+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


++++++++++++++++++H++++++++4++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


POLLUTANT LOADING; DIVIDECK BASIN 
EVENTS of: JANUARY 
RPC CONCENTRATION DATA 


LAND use: rural- 0.921 LAC- 0.011 LAR* O.O&I MAC- 0.001 HAR- O.OOf CPD* 0.00 


POLLUTANT! 

LOAD-KOS. 


SUSPENDED 

SOLIDS 


0.000 


BOD 


0.000 


COD 


0.000 


TOTAL KJELDAHL 
NITROGEN 


0.000 


TOTAL 

PHOSPOROUS 


0.000 


COG CONCENTRATION DATA 


LAND use; shop CTR- 0.00 1 STRP COM- 0.01 » HR RES- 
LLSFR- 0.06»CON» O.OOf FOR= 0.82f CT 


O.OOf TH/GA' 
AG- O.OOf MT 


AG> 


,00f HDSFR- 
O.OOf PAST 


pollutant; 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


0.000 


0.00 


TOTAL KJELDAHL 
NITROGEN 


0.000 


TOTAL 

NITROGEN 


0.000 


TOTAL 

PHOSPOROUS 


0.000 


+++4++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4444444444444444444444444444444444444444444444444+++++++++++++++++++++++++++++++++++++++++++++++++++ 


0.00 
' 0.10 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++t++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


++++++++f++++++++++++1-++++ttt +++++++++++++++++++ +++++t ++++++++++++++ t+++t+++++++++ ++++++++++++++++++ 


4 

4 

4 

4- 

TIME 

DISCHARGE 

RUNOFP SUMMARY FOR DIVIDECK BASIN 
AUG. EUENT FOR: FEBRUARY 
AUG. NO, OF storms: 7.63 

RAIN RATE CUMULATIUE 

CUMULATIUE 

+ 

+ 

+ 

+ 

+ 

+ 

4- 

4- 

(SECS) 

(M3/S) 


(M/S) 

DISCHARGE(M3> 

RAINFALL(M3) 

4 

4 

+ 

4- 

900. 

0.00000000 


0.0000000420 

0.0000 

149.3732 

4 

4 

+ 

1800. 

0.00000000 


0,0000000420 

0.0000 

290.7465 

4 

+ 

2700. 

0,00000000 


0.0000001250 

0.0000 

743.3097 

4 

+ 

3600. 

0.00000000 


0.0000001250 

0,0000 

1187.8729 

4 

4- 

4500. 

0,00000000 


0.0000002220 

0.0000 

1977.4172 

4 

4- 

5400. 

0.00000000 


0.0000002220 

0.0000 

2766.9614 

4 

4- 

6300. 

0.00000000 


0.0000003060 

0.0000 

3855.2522 

4 

4- 

7200. 

0.00000000 


0.0000003060 

0,0000 

4943.5430 

4 

4- 

8100. 

0.00000000 


0.0000003340 

0.0000 

6131.4160 

4 

4- 

9000. 

0.00000000 


0.0000003340 

0.0000 

7319.2886 

4 

4- 

9900. 

0.00000000 


0.0000004720 

0.0000 

8997.9590 

4 

4- 

10800. 

0.00000000 


0.0000004720 

0.0000 

10676.6299 

4 

+ 

11700. 

0,00000000 


0.0000005560 

0,0000 

12654.0469 

4 

4- 

12600, 

0.00000000 


0.00000055ib0 

0.0000 

14631.4639 

4 

4- 

13500. 

0,00000000 


0.0000006670 

0.0000 

17003.6523 

4 

4- 

14400. 

0.00000000 


0.0000006670 

0.0000 

19375.8418 

4 

4- 

15300, 

0.00000000 


0.0000006110 

0.0000 

21548.8672 

4 

4- 

16200. 

0.00000000 


0.0000006110 

0.0000 

23721.8926 

4 

4 

17100. 

0.00000000 


0.0000005280 

0.0000 

25599,7285 

4 

4 

18000. 

0.00000000 


0.0000005280 

0.0000 

27477.5625 

4 

4 

18900. 

0.00001425 


0.0000004440 

0.0128 

29056.6504 

4 

4 

19800. 

0.00033550 


0,0000004440 

0.3148 

30635.7402 

4 

4 

20700. 

0,00169698 


0,0000003610 

1.8421 

31919,6387 

4 

4 

21600, 

0.00542458 


0.0000003610 

6.7242 

33203,5391 

4 

4 

22500. 

0,01116251 


0.0000002780 

16.7705 

34192.2461 

4 

4 

23400. 

0.01898889 


0.0000002780 

33.8605 

35180.9570 

4 

4 

24300. 

0.02618044 


0,0000001940 

57.4229 

35870.9180 

4 

4 

25200. 

0.03221604 


0,0000001940 

86.4173 

36560.8828 

4 

4 

26100. 

0.03465437 


0.0000001110 

117,6063 

36955.6562 

4 

4 

27000, 

0.03447625 


0.0000001110 

148.6349 

37350.4258 

4 

4 

27900. 

0.03222093 


0.0000000420 

177.6337 

37499.8008 

4 

4 

28000. 

0.02927647 


0,0000000420 

203.9025 

37649.1719 

4 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++4+++++++++++++++++++++++++++++++++4+++++++++++++++ 
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+++++++++++++++++++++++++++++++++++++++4+++4++++++++++++++++++++4+++++4+++++++++4+++++++++++++++++++ 


4444444444444444444444444444444444444444444444444444444444444444444444+44444444444444444444444444444 


POLLUTANT LOADING: DIVIDECK BASIN 
EVENTS of; FEBRUARY 
RPC CONCENTRATION DATA 


LAND use: RURAL=> 0.92( LAC=c 0.01$ LAR» 0.06) HAC> 0.00) HAR=> 0.00) CBD> 0.00 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


274.036 


BOD 


10.978 


COD 


63.728 


total KJELDAHL 
NITROGEN 


1.271 


TOTAL 

PHOSPOROUS 


0.647 


COG CONCENTRATION DATA 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


99.254 


COD 


190.49 


TOTAL KJELDAHL 
NITROGEN 


2.704 


TOTAL TOTAL 

NITROGEN PHOSPOROUS 


2.855 


0.356 


4444+4444444444444444444+4+4+++++++++++++++++++++44+++444444444+4+4+44444444+4444444+444+44444444444 
444444444444444444444444444+444+++44++4++444444444++4+44444444444+4444444444444444444444444444444444 
4 


4 
4 
4 
4 
4 

LAND use: SHOP CTR- 0.00) STRP COM= 0.01) HR RES- 0,00) TH/GA- 0.00) MDSFR- O'.OO 4 
LLSFR- 0.06) CON- 0.00) FOR- 0.82) CT AG- 0.00) MT AG- 0.00) PAST- 0.10 4 

4 
4 
4 
4 
4 
4 
4 
4 
4 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

44444444444444444444444444444444++++44+4+4++44+++444444+444444444444444*44444444444444444444444444444 
44+444444444444+4444444444444444+++++++44+4+44++4444+444+4444444444+44444++4++444444444 4444444444444 
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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


T T X T T 

4 

4 

4 

4 

4 

4 

4 

4 

TIME 

DISCHARGE 

RUNOFF 

AVG 

AUG. 

SUMMARY FOR DIUIDECK BASIN 
. EUENT for: march 
NO. OF storms: 8.78 

RAIN RATE CUMULATIUE 

CUMULATIUE 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4 

4 

(SECS) 

<M3/S) 


(M/S) 

DISCHARGE(M3) 

RAINFALL(M3) 

+ 

+ 

4 

4 

900. 

0.00000000 


0.0000000520 

0.0000 

184.9383 

+ 

+ 

4 

1800. 

0.00000000 


0.0000000520 

0.0000 

369.8766 

+ 

4 

2700. 

0.00000000 


0.0000001570 

0.0000 

928.2480 

+ 

4 

3600. 

0.00000000 


0,0000001570 

0.0000 

1486.6194 

+ 

4 

4500. 

0.00000000 


0.0000002620 

0,0000 

2418.4238 

+ 

4 

5400. 

0.00000000 


0.0000002620 

0,0000 

3350.2283 

4 

4 

6300. 

0.00000000 


0.0000003670 

0.0000 

4655.4658 

4 

4 

7200. 

0.00000000 


0.0000003670 

0.0000 

5960.7036 

4 

4 

0100. 

0.00000000 


0.0000004710 

0.0000 

7635.8179 

4 

4 

9000. 

0.00000000 


0.0000004710 

0.0000 

9310.9316 

4 

4 

9900. 

0.00000000 


0.0000005760 

0.0000 

11359.4785 

4 

4 

10800. 

0.00000000 


0.0000005760 

0.0000 

13408.0264 

4 

4 

11700. 

0.00000000 


0.0000006390 

0.0000 

15680.6338 

4 

4 

12600. 

0.00000000 


0.0000006390 

0.0000 

17953.2402 

4 

4 

13500. 

0.00000000 


0.0000006210 

0.0000 

20161.8301 

4 

4 

14400. 

0.00000000 


0.0000006210 

0.0000 

22370.4199 

4 

4 

15300. 

0.00000000 


0.0000005460 

0.0000 

24312.2715 

4 

4 

16200. 

0.00000701 


0,0000005460 

0.0063 

26254.1250 

4 

4 

17100. 

0.00030524 


0,0000004710 

0.2810 

27929.2383 

4 

4 

IBOOO. 

0.00271932 


0.0000004710 

2.7284 

29604,3535 

4 

4 

10900. 

0.01004917 


0.0000004040 

11.7727 

31041.1816 

4 

4 

19800. 

0.02589324 


0.0000004040 

35.0766 

32478,0090 

4 

4 

20700. 

0.04984923 


0.0000003370 

79.9409 

33676,5508 

4 

4 

21600. 

0.08218576 


0.0000003370 

153.9081 

34875.0937 

4 

4 

22500. 

0.11506171 


0.0000002620 

257.4636 

35806.8984 

4 

4 

23400. 

0.14304954 


0.0000002620 

386.2082- 

36738,7031 

4 

4 

24300. 

0.15575710 


0.0000001870 

526.3896 

37403.7695 

4 

4 

25200. 

0.15440488 


0.0000001870 

665.3541 

38068.8359 

4 

4 

26100. 

0.13917698 


0.0000001120 

790.6133 

38467.1641 

4 

4 

27000. 

0.11739077 


0.0000001120 

896.2650 

38865.4922 

4 

4 

27900. 

0.09321304 


0.0000000370 

980.1567 

38997,0820 

4 

4 

28800. 

0.07302280 


0.0000000370 

1045.8773 

39128.6719 

4 


+ + 

+ + 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

+++++++++++++++++++++++++4+++++++++++++++++++++++++++++++++++++++++++++++++++++++4++++++++++++++++++ 
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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


+++++++++++++++++++++++++++++++++++++++++++++++4++++++++++++++++++++++++++++++++++++++++++++++++++++ 


POLLUTANT LOADING: DIUIDECK BASIN 
EVENTS of: march 
RPC CONCENTRATION DATA 


LAND use: RURAL* 0.92r LAC= O.OH LAR* 0.06» HAC- O.OOf HAR= O.OOi CBD* 0.00 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


1616.834 


BOD 


COD 


64.768 376.001 


TOTAL KJELDAHL 
NITROGEN 


7.S07 


TOTAL 

PHOSPDROUS 


3.820 


COG CONCENTRATION DATA 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


585.606 


4444444444444444444444444444444+444444444444+444444444444444444444+44+444444444444444444444444444444 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++44+++++++++++++++++++++++++++++++++++++++++++ 


+ 
+ 
+ 
+ 
+ 

LAND use: shop CTR= O.OOf STRP COM= 0.01 r HR RES* O.OO: TH/GA= 0.00» MDSFR= 0.00 + 

LLSFR= 0.06: CON= O.OOf FOR= O.B2» CT AQ= O.OOf MT AG= O.OOf PAST* 0.10 + 

+ 

4 

COD TOTAL KJELDAHL TOTAL TOTAL 4 

NITROGEN NITROGEN PHOSPOROUS 4 

4 
4 

1123.89 15.951 16.846 2.102 4 

4 
4 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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+4 ++++++++++++++++++++++++++++++++++++++++++++++++4+++++++++++++++++++++++++++++++++++++++++++++++++ 
4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


4 

4 

4 

4 

4 

4 

TIME 

DISCHARGE 

RUNOFF 

AUG 

AUG. 

SUMMARY FOR DIUIDECK BASIN 
. EVENT for; APRIL 
NO. OF storms: 10.24 

RAIN RATE CUMULATIVE 

CUMULATIVE 

+ 

+ 

4 

4 

4 

4 

4 

4 

(SECS) 

(M3/S) 


(M/S) 

DISCHAR6E(M3) 

RAINFALL(M3) 

4 

4 

4 

4 

900. 

0.00000000 


0.0000000510 

0.0000 

181.3818 

4 

4 

4 

1800. 

0.00000000 


0.0000000510 

0.0000 

362.7636 

4 

4 

2700, 

0.00000000 


0,0000001660 

0.0000 

953.1436 

4 

4 

3600. 

0.00000000 


0.0000001660 

0.0000 

1543.5234 

4 

4 

4500, 

0.00000000 


0.0000002750 

0.0000 

2521.5625 

4 

4 

5400. 

0.00000000 


0.0000002750 

0.0000 

3499.6016 

4 

4 

6300. 

0.00000000 


0.0000003900 

0.0000 

4886.6387 

4 

4 

7200. 

0.00000000 


0.0000003900 

0.0000 

6273.6763 

4 

4 

8100. 

0.00000000 


0.0000004980 

0.0000 

8044.8164 

4 

4 

9000. 

0.00000000 


0.0000004980 

0,0000 

9815.9561 

4 

4 

9900. 

0.00000000 


0.0000006070 

0.0000 

11974.7549 

4 

4 

10800. 

0.00000000 


0,0000006070 

0.0000 

14133.5537 

4 

4 

11700. 

0.00000000 


0.0000007220 

0,0000 

16701.3516 

4 

4 

12600. 

0.00000000 


0.0000007220 

0.0000 

19269.1484 

4 

4 

13500. 

0.00000000 


0.0000005240 

0.0000 

21132,7578 

4 

4 

14400. 

0.00000000 


0,0000005240 

0.0000 

22996.3672 

4 

4 

15300. 

0.00000000 


0,0000004280 

0.0000 

24518.5508 

4 

4 

16200. 

0,00000000 


0.0000004280 

0.0000 

26040.7363 

4 

4 

17100. 

0.00000000 


0.0000003320 

0,0000 

27221.4961 

4 

4 

18000. 

0.00000000 


0.0000003320 

0.0000 

28402.2559 

4 

4 

18900. 

0.00000005 


0.0000002360 

0.0000 

29241.5918 

4 

4 

19800. 

0.00000057 


0.0000002360 

0.0006 

30080.9258 

4 

4 

20700, 

0.00000057 


0.0000001470 

0.0011 

30603.7324 

4 

4 

21600. 

0.00000057 


0.0000001470 

0.0016 

31126.5391 

4 

4 

22500. 

0,00000057 


0.0000000510 

0.0021 

31307.9199 

4 

4 

23400. 

0.00000057 


0.0000000510 

0 ^ 0P26 

31489.3027 

4 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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444444444444444444444444444444444444444+444444444444444444444444444444444444444444444444444444444444 


POLLUTANT LOADING: DIVIDECK BASIN 
EVENTS of: APRIL 
RPC CONCENTRATION DATA 


LAND use; RURAL= 0.92? LAC= O.Olr LAR= O.OAf HAC= 0.005 HAR» 0.005 CBD= 0.00 


pollutant; 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


0.005 


BOD 


0.000 


COD 


0.001 


TOTAL KJELDAHL 
NITROGEN 


0.000 


TOTAL 

FHOSPOROUS 


0.000 


COG CONCENTRATION DATA 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


0.002 


COD 


0.00 


TOTAL KJELDAHL 
NITROGEN 


0.000 


TOTAL TOTAL 

NITROGEN PHOSPOROUS 


0.000 


0.000 


44444444+44+44444444444444++444444444+44+44444444444444+44444+44444444444444444444444444444444444444 
44444444444+44+4444444+44444444+4444444444444++444+4444++44+4+444444+4444444444++4444444+44444444444 
4 


4 
4 
4 
4 
4 

LAND use: SHOP CTR= 0.005 STRP COM= 0.015 HR RES- 0.005 TH/QA- 0.005 MDSFR- 0.00 4 

LLSFR=' 0.065 CON= 0.005 FOR= 0.825 CT AG- 0.005 MT AG- 0.005 PAST- 0.10 4 

4 
4 
4 
4 
4 
4 
4 
4 
4 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4444444444444444444444444444444444444444444+444444444444444444444S- 4+44444444444444444444444444444444 
4444444+444444444444444444444444444+44+444444444444444+44444+4++44444444444+4444+444444444+4+4444444 
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++++++++++++++++++++++++++++++++++++++++++++++ 4 -+++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


RUNOFF SUMMARY FOR DIVIDECK BASIN 
AUG. EVENT FOR: MAY 
AVG. NO. OF storms: 10.52 


+ 

TIME 

PISCHARGE 

RAIN RATE 

CUMULATIVE 

CUMULATIVE 

+ 

+ 

+ 

<SECS) 

<M3/S) 

(M/S) 

DISCHARGE(M3) 

RAINFALL(M3) 

+ 

+ 

+ 

+ 

900. 

0,00000000 

0.0000000940 

0.0000 

334.3115 

+ 

+ 

+ 

IBOO, 

0.00000000 

0.0000000940 

0.0000 

668.6230 

+ 

+ 

2700. 

0.00000000 

0.0000002810 

0.0000 

1668.0012 

+ 

+ 

3600. 

0.00000000 

0.0000002810 

0.0000 

2667.3794 

+ 

+ 

4500. 

0.00000000 

0.0000004690 

0.0000 

4335.3809 

+ 

+ 

5400. 

0.00000000 

0.0000004690 

0.0000 

6003.3816 

+ 

+ 

6300. 

0.00000000 

0.0000006570 

0.0000 

8340.0059 

+ 

+ 

7200. 

0.00000000 

0.0000006570 

0,0000 

10676.6299 

+ 

+ 

8100. 

0,00000000 

0.0000008310 

0.0000 

13632.0859 

+ 

+ 

9000. 

0.00000000 

0.0000008310 

0.0000 

16587.5430 

+ 

+ 

9900. 

0.00000000 

0.0000010000 

0.0000 

20144.0468 

+ 

+ 

loeoo. 

0.00000000 

0.0000010000 

0,0000 

23700.5547 

+ 

+ 

11700. 

0.00000000 

0.0000007100 

0.0000 

26225.6738 

+ 

+ 

12600. 

0.00000000 

0,0000007100 

0.0000 

28750.7930 

+ 

+ 

13500. 

0.00019675 

0.0000004960 

0.1771 

30514.8203 

+ 

+ 

14400. 

0.00307646 

0.0000004960 

2.9459 

32278.8457 

+ 

+ 

15300. 

0.00918785 

0.0000002950 

11.2150 

33328.0156 

+ 

+ 

16200. 

0.01883091 

0.0000002950 

26.1628 

34377.1836 

+ 

+ 

17100. 

0.02505752 

0.0000000940 

50.7145 

34711.4961 

+ 

+ 

IBOOO. 

0.02774033 

0.0000000940 

75.6808 

35045.8086 

+ 


+++++++++++++++++++++++++++++++++++++t++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


POLLUTANT LOADING.* DIUIDECK BASIN 
EUENTS of: may 
RPC CONCENTRATION DATA 


LAND use: RURAL= 0.92f LAC= O.OIJ LAR= 0.06» HAC= 0.00» HAR= O.OOS CBD= 0.00 


pollutant: 

LOAD-KQS. 


SUSPENDED 

SOLIDS 


1A0.182 


BOD 


5.615 


COD 


32.600 


TOTAL KJELDAHL 
NITROGEN 


0.651 


TOTAL 

PHOSPOROUS 


0.331 


pollutant: 

LOAD-KGS, 


SUSPENDED 

SOLIDS 


50.773 


COD 


97.44 


TOTAL KJELDAHL 
NITROGEN 


1.383 


TOTAL 

NITROGEN 


1.461 


TOTAL 

PHOSPOROUS 


0.182 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
4 


COG CONCENTRATION DATA 


LAND use: SHOP CTR=* O.OOJ STRP COM- 0.011 HR RES- 0.001 TH/GA- 0.001 MDSFR- 0.00 

LLSFR= 0.06» CON- 0.00» FOR= 0.821 CT AG- O.OOf MT AG= 0.00» PAST- 0.10 


4444444444444444444444444444444444444444444444444444444444444444444444+44-444444444444444444444444444 

44444444444444444444444444+44444444444444444444+4444444444444444444444+44444444444444+4444444444+444 


.-17 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


RUNOFF SUMMARY FOR DIVIDECK BASIN 
AUG. EUENT FOR: JUNE 
AUG. NO. OF storms: 9.A5 


+ 

TIME 

DISCHARGE 

RAIN RATE 

CUMULATIUE 

CUMULATIUE 

+ 

+ 

+ 

(SECS) 

(M3/S) 

(M/S) 

DISCHARGE(M3) 

RAINFALL(M3) 

+ 

+ 

+ 

+ 

900. 

0.00000000 

0.0000001100 

0.0000 

391 .2156 

+ 

+ 

+ 

IBOO. 

0.00000000 

0.0000003290 

0.0000 

1561.3060 

+ 

+ 

2700. 

0.00000000 

0.0000005490 

0.0000 

3513.8276 

+ 

+ 

3600. 

0.00000000 

0.0000007680 

0.0000 

6245.2241 

+ 

+ 

4500. 

0.00000000 

0.0000009880 

0.0000 

9759.0518 

+ 

+ 

5400. 

0.00000000 

0.0000012070 

0.0000 

14051.7539 

+ 

+ 

6300. 

0.00000000 

0.0000013200 

0.0000 

18746.3418 

+ 

+ 

7200. 

0.00000000 

0.0000013200 

0.0000 

23440.9297 

+ 

+ 

BIOO. 

0.00000000 

0.0000011440 

0.0000 

27509.5723 

+ 

+ 

9000. 

0.00000000 

0.0000009560 

0.0000 

30909.5918 

+ 

+ 

9900. 

0.00089749 

0.0000007370 

0.8077 

33530.7383 

+ 

+ 

10800. 

0.00787075 

0.0000005170 

7.8914 

35369.4531 

+ 

+ 

11700. 

0.02187117 

0.0000003290 

27.5755 

36539.5430 

+ 

+ 

12600. 

0.033271B7 

0.0000001100 

57.5201 

36930.7578 

+ 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


POLLUTANT loading: divideck dasin 
EVENTS of; JUNE 
RPC CONCENTRATION DATA 


LAND use; rural* 0.92t LAC= O.Oi; LAR= O.O&J HAC= O.OOf HAR= 0.00) CBD= 0.00 


pollutant; 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


97.732 


BOD 


COD 


22,728 


TOTAL KJELDAHL 
NITROGEN 


TOTAL 

PHOSPOROUS 


0.231 


COG CONCENTRATION DATA 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


35,398 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
+++++++++++++++++++++++++++++++++++4++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
4 


4 
4 
4 
4 
4 

LAND use; SHOP CTR* 0.00) STRP COM= 0.01) HR RES* 0.00) TH/GA* 0,00) MDSFR* 0.00 4 

LLSFR= 0.06) CON= 0.00) FOR* 0,82) CT AG= 0.00) MT AG= 0.00) PAST* 0.10 4 

4 

4 

COD TOTAL KJELDAHL TOTAL TOTAL 4 

NITROGEN nitrogen PHOSPOROUS 4 

4 
4 

67.94 0.964 I.OIB 0.127 4 

4 
4 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
+ 

444444444444444444444444444444444444444444444444444444444444444444+4'44444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


H 

VO 
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444444444444444444444444444444444+44444444444444++++++++++++444444444444444444444444444444444444++++ 


RUNOFF SUMMARY 
AUG. EUENT 
AUG. NO. OF 


FOR DIUIDECK 
FOR! JULY 
storms: 9 


11 


4 

TIME 

DISCHARGE 

RAIN RATE 

CUMULATIUE 

CUMULATIVE 

+ 

4 

4 

(SECS) 

(M3/S) 

(M/S) 

DISCHARGE(M3> 

RAINFALL(M3) 

+ 

+ 

4 

4 

900. 

0.00000000 

0.0000001230 

0.0000 

437.4502 

+ 

4 

1000. 

0.00000000 

0.0000003540 

0.0000 

1696.4532 

+ 

4 

2700. 

0.00000000 

0.0000005850 

0.0000 

3777.0093 

+ 

4 

3600. 

0.00000000 

0.0000008160 

0.0000 

6679.1182 

+ 

4 

4500. 

0.00000000 

0.0000010480 

0.0000 

10406.3359 

+ 

4 

5400. 

0.00000000 

0.0000012790 

0.0000 

14955.1064 

+ 

4 

6300. 

0.00000000 

0.0000015100 

0.0000 

20325,4297 

+ 

4 

7200. 

0.00000000 

0.0000016940 

0.0000 

26350.1504 

+ 

4 

8100. 

0.00036758 

0.0000013250 

0.3308 

31062.5195 

+ 

4 

9000. 

0.02616293 

0.0000011050 

23.8775 

34992.4570 

+ 

4 

9900. 

0.13193822 

0.0000008430 

142.6219 

37990.5898 

+ 

4 

10800. 

0.31343022 

0.0000005980 

424.7090 

40117.3789 

+ 

4 

11700. 

0.49779BB3 

0.0000003670 

872.7280 

41422.6172 

+ 

4 

12600. 

0.57992494 

0.0000001230 

1394.6605 

41860.0664 

+ 


4444++++++4+ h+4444444++44444444444444+++ + + ++++ + + 444444 +++++++ + + ++4+++4444444444+444444444444 + 4444444 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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POLLUTANT LOAPINGI DIVIDECK BASIN 
EVENTS of: JULY 
RPC CONCENTRATION DATA 


LAND use: RURAL= 0.92t LAC= 0.015 LAR= 0.065 MAC- 0.005 HAR- 0.005 CBD- 0.00 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


2237.058 


BOD 


89.613 


COD 


520.236 


TOTAL KJELDAHL 
NITROGEN 


10.387 


TOTAL 

PHOSPOROUS 


5. 285 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


010.246 


COD 


1555.02 


TOTAL KJELDAHL 
NITROGEN 


22.070 


TOTAL 

NITROGEN 


23.308 


TOTAL 

PHOSPOROUS 


2.909 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


COG CONCENTRATION DATA 


LAND use: SHOP CTR«= 0.005 STRP COM- 0.015 HR RES- 0.005 TH/GA- 0.005 MDSFR- 0.00 

LLSFR- 0.065 CON= 0.005 FOR- 0.825 CT AO- 0.005 MT AG- 0.005 PAST- 0.10 


+++++++++++4+++44444444444444444444444444444444444444444444444444V4444444444444444444444444444444444 

+++++44+++++++++++++++++++4++++++4444444444444444444444444444444444444444444444444444444444444444444 
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4 4 
4 4 
4 RUNOFF SUMMARY FOR DIVIDECK BASIN 4 
4 AVG. EVENT FOR: AUGUST 4 
4 AVG. NO. OF storms; 9.46 4 
4 4 
4 4 


4 

TIME 

DISCHARGE 

RAIN RATE 

CUMULATIVE 

CUMULATIVE 

+ 

4 

4 

(SECS) 

<M3/S) 

(M/S) 

DISCHARGE(M3) 

RAINFALL<M3) 

+ 

+ 

4 

4 

900. 

0.00000000 

0.0000001080 

0.0000 

364.1026 

+ 

+ 

4 

1600. 

0.00000000 

0.0000003230 

0.0000 

1532.8540 

+ 

4 

2700. 

0.00000000 

0.0000005390 

0.0000 

3449.8105 

+ 

4 

3600. 

0.00000000 

0.0000007700 

0,0000 

6188.3198 

+ 

4 

4500. 

0.00000000 

0.0000010010 

0.0000 

9748.3818 

+ 

4 

5400. 

0.00000000 

0.0000012170 

0.0000 

14076.6494 

+ 

4 

6300. 

0.00000000 

0.0000014320 

0.0000 

19169.5664 

+ 

4 

7200. 

0.00000000 

0.0000015950 

0.0000 

24842.1934 

+ 

4 

8100. 

0.00000000 

0.0000015950 

0.0000 

30514.8203 

+ 

4 

9000. 

0.00278445 

0.0000014630 

2,5060 

35717,9883 

+ 

4 

9900. 

0.06634991 

0.0000012320 

62.2209 

40099.6055 

+ 

4 

10600. 

0.27713221 

0.0000010010 

311.6399 

43659.6680 

+ 

4 

11700. 

0.61399269 

0.0000007850 

864.2334 

46451.5234 

+ 

4 

12600. 

0.92908120 

0.0000005540 

1700.4065 

48421.8281 

4 

4 

13500. 

1.06917965 

0.0000003390 

2662.6682 

49627.4844 

+ 

4 

14400. 

1.00847340 

0.0000001080 

3570.2942 

50011.5859 

+ 


4 


4 


4 


4 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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44444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 ^^ 


POLLUTANT LOADING! DIVIDECK BASIN 

EVENTS of: august 

RPC CONCENTRATION DATA 


LAND use: RURALb 0.92i LAC= 0.01 r LAR» 0.06r HAC>= O.OOt HAR> O.OOf CBD> 0.00 

BOB 


pollutant: 

LOAD-KQS. 


SUSPENDED 

SOLIDS 


5946,829 


COD 


238,222 1382.957 


TOTAL KJELDAHL 
NITROGEN 


27.612 


TOTAL 

PHOSPOROUS 


14,049 


COG CONCENTRATION DATA 


POLLUTANT ! 
LOAD-KGS. 


SUSPENDED 

SOLIDS 


COD 


2153.898 4133.74 


TOTAL KJELDAHL 
NITROGEN 


58.668 


TOTAL 

NITROGEN 


61.960 


TOTAL 

PHOSPOROUS 


7.732 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


4 
4 
4 
4 
4 

LAND USE! SHOP CTR« 0.00! STRP COM= 0.011 HR RES= 0.00! TH/GA- 0.00! MDSFR= 0,00 4 

LLSFR= 0.061 CON* 0.00! FOR= 0.82! CT AG* 0.001 MT AG* 0.00! PAST* 0.10 4 

4 
4 
4 
4 
4 
4 
4 
4 
4 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
4 


RUNOFF SUMMARY 
AUG. EUENT 
AUQ. NO. OF 


FOR DIUIDECK BASIN 

for: septembr 
storms; 7.29 


4 

TIME 

DISCHARGE 

RAIN RATE 

CUMULATIUE 

CUMULATIVE 

+ 

4 

4 

(SECS) 

(M3/S) 

(M/S) 

niSCHARGE(M3) 

RAINFALL(M3) 

+ 

+ 

4 

4 

900. 

0.00000000 

0.0000001040 

0.0000 

369.8766 

+ 

+ 

4 

1800. 

0.00000000 

0.0000001040 

0.0000 

739.7532 

+ 

4 

2700. 

0.00000000 

0.0000003120 

0.0000 

1849.3829 

+ 

4 

3600. 

0.00000000 

0.0000003120 

0.0000 

2959.0127 

+ 

4 

4500. 

0.00000000 

0,0000005050 

0.0000 

4755.0479 

+ 

4 

5400. 

0.00000000 

0.0000005050 

0.0000 

6551 . 0835 

+ 

4 

6300. 

0.00000000 

0,0000007130 

0.0000 

9086.8721 

+ 

4 

7200. 

0.00000000 

0.0000007130 

0.0000 

11622.6611 

4 

4 

8100. 

0.00000000 

0.0000009210 

0.0000 

14898.2031 

4 

4 

9000. 

0.00000000 

0,0000009210 

0.0000 

18173.7441 

4 

4 

9900. 

0.00000000 

0.0000009900 

0.0000 

21694.6855 

4 

4 

10800. 

0.00000000 

0.0000009900 

0.0000 

25215.6250 

4 

4 

11700. 

0.00000000 

0.0000008320 

0.0000 

28174.6387 

4 

4 

12600. 

0.00000000 

0,0000008320 

0.0000 

31133.6504 

4 

4 

13500. 

0.00067487 

0.0000006530 

0.6074 

33456.0469 

4 

4 

14400. 

0.01051014 

0.0000006530 

10.0665 

35778.4453 

4 

4 

15300. 

0^03753251 

0.0000004600 

43.8458 

37414.4375 

4 

4 

16200. 

0.08861833 

0.0000004600 

123.6023 

39050.4297 

4 

4 

17100. 

0.14782076 

0.0000002820 

256.6410 

40053.3633 

4 

4 

18000. 

0,20524813 

0.0000002820 

441.3643 

41056.3008 

4 

4 

18900. 

0.22742464 

0.0000000890 

646.0464 

41372.8281 

4 

4 

19800. 

0.21599340 

0.0000000890 

840.4405 

41689.3594 

4 


444444444-444444444444444444444444444444444'44444444444444444444444444444-44444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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444444444444444444444444444+444++444++44+44444+44++44444++4+44+++4444++++++44444+4444+444444+444444H 


POLLUTANT LOADING: DIVIDECK BASIN 
EVENTS of: SEFTEMBR 
RPC CONCENTRATION DATA 


LAND use: rural- 0.921 LAC- 0.01 f LAR= 0.06» MAC- 0.001 HAR- O.OOJ CBD= 0.00 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


1078.760 


BOD 


43.214 


COD 


250.870 


TOTAL KJELDAHL 
NITROGEN 


5.009 


TOTAL 

PHOSFOROUS 


2.549 


COG CONCENTRATION DATA 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


390.719 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


4 
4 
4 
4 
4 

LAND use: shop CTR= O.OOS STRP C0M= 0.01 » HR RES= O.OOt TH/GA- 0.001 MDSFR= 0.00 4 

LLSFR- 0.061 CON- 0.00! FOR- 0.82! CT AG= 0.00! MT AG= 0.001 FAST- 0.10 4 

4 
4 

COD TOTAL KJELDAHL TOTAL TOTAL 4 

NITROGEN NITROGEN PHDSPOROUS 4 

4 
4 

749.86 10.643 11.240 1.403 4 

4 

4 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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4444444444+4444444444+44444+444444444+444444+4444+444444444444+4444+44444444444444444444444^44444444 


4 

4 

4 

4 

4 

4 

TIME 

DISCHARGE 

RUNOFF 

AVG 

AVG. 

SUMMARY FOR DI VI DECK BASIN 
. EVENT for: OCTOBER 
NO. OF storms: 6,60 

RAIN RATE CUMULATIVE 

CUMULATIVE 

4 

4 

4 

4 

4 

4 

4 

4 

(SECS) 

(M3/S) 


(M/S) 

DISCHARGE(M3) 

RAINFALL(M3) 

■f 

4 

4 

4 

900. 

0.00000000 


0.0000000620 

0.0000 

220.5033 

4 

4 

4 

1800. 

0.00000000 


0.0000000620 

0.0000 

441.0067 

4 

4 

2700. 

0.00000000 


0.0000001860 

0.0000 

1102.5167 

4 

4 

3600. 

0.00000000 


0.0000001860 

0.0000 

1764,0267 

4 

4 

4500. 

0.00000000 


0.0000003090 

0.0000 

2862.9871 

4 

4 

5400. 

0.00000000 


0.0000003090 

0.0000 

3961.9473 

4 

4 

6300. 

0.00000000 


0.0000004330 

0.0000 

5501.9141 

4 

4 

7200. 

0.00000000 


0.0000004330 

0.0000 

7041 .8813 

4 

4 

8100. 

0.00000000 


0.0000005570 

0,0000 

9022.8555 

4 

4 

9000. 

0.00000000 


0.0000005570 

0.0000 

11003.8291 

4 

4 

9900. 

0.00000000 


0.0000006810 

0.0000 

13425.8096 

4 

4 

10800. 

0.00000000 


0.0000006810 

0.0000 

15847.7900 

4 

4 

11700. 

0.00000000 


O.OOOOOOBBOO 

0.0000 

18977,5156 

4 

4 

12600. 

0,00000000 


0.0000008800 

0.0000 

22107.2402 

4 

4 

13500. 

0.00000000 


0.0000007300 

0.0000 

24703,4902 

4 

4 

14400. 

0.00000000 


0.0000007300 

0.0000 

27299.7383 

4 

4 

15300. 

0.00000000 


0.0000005940 

0.0000 

29412,3027 

4 

4 

16200. 

0.00000000 


0.0000005940 

0.0000 

31524.8672 

4 

4 

17100. 

0.00000000 


0.0000005070 

0,0000 

33328,0156 

4 

4 

18000, 

0,00014296 


0.0000005070 

0,1287 

35131 . 1641 

4 

4 

18900. 

0.00156652 


0.0000004080 

1.5385 

36582.2187 

4 

4 

19800. 

0.00692867 


0.0000004080 

7.7743 

38033.2734 

4 

4 

20700. 

0.01691719 


0.0000003220 

22.9998 

39178.4687 

4 

4 

21600. 

0.03244998 


0.0000003220 

52.2048 

40323.6641 

4 

4 

22500, 

0.04948571 


0.0000002350 

96.7419 

41159,4414 

4 

4 

23400. 

0.06656066 


0.0000002350 

156.6465 

41995.2227 

4 

4 

24300, 

0.07638235 


0.0000001360 

225.3906 

42478.9062 

4 

4 

25200. 

0.07803639 


0,0000001360 

295,6234 

42962.5937 

4 

4 

26100. 

0.07063751 


0.0000000500 

359.1971 

43140.4180 

4 

4 

27000. 

0.06032138 


0.0000000500 

413.4864 

43318.2461 

4 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


+ + 

+ + 

+ POLLUTANT loading; niVIDECK BASIN + 

+ EVENTS OF! OCTOBER + 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4 

4 

4 

4 

4 


RPC CONCENTRATION DATA 

LAND USE! RURAL= 0.92J LAC= O.Olf LAR= 0.06? HAC= 0.00? HAR= 0.00? CBD= 0.00 

POLLUTANT! SUSPENDED BOD COD TOTAL KJELDAHL TOTAL 

LOAD-KGS. SOLIDS NITROGEN PHOSPOROUS 

400.502 19.240 111.743 2.231 1.135 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


4444+44444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


++ 44 +++ 4 +++++++++++++++++++++++++++++++++++ 44 ++++++++++++++++++++++++++ 4 ++++++++++++ 4 + 44+++++++++++4 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


COG CONCENTRATION DATA 

LAND USE! SHOP CTR= 0.00? STRP COM= 0.01? HR RES= 0.00? TH/GA= 0.00? HDSFR* 0.00 

LLSFR- 0.06? CON= 0.00? FOR= 0.82? CT AG* 0.00? HT AG= 0.00? PAST* 0,10 

pollutant: suspended cod total KJELDAHL TOTAL TOTAL 

LOAD-KGS, SOLIDS NITROGEN NITROGEN PHOSPOROUS 

174.034 334.01 4,740 5.006 0.625 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


4 4 
44 +++ 444 +++ 4 ++++ 44 +++++++++++++++++ 44 + 4 +++++++ 44 ++++++ 44 ++++++++++ 4444 ++++++ 444 ++ 444 ++++++ 4 + 4 +++ 4++4 
++++ 44 +++++++ 4 ++ 4 +++++++++++++++++++++++++++++++++++ 4 ++ 4 +++++++ 444 ++++ 4 +++ 4444 + 44 +++++ 4 ++++ 4 ++++++++ 
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+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ++ + + + + + + 4 + + + + + + + + + + + + 4 + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 4 + -t ^ 


■» 

4 

4 

4 

4 

4 

TIME 

DISCHARGE 

RUNOFF 

AOG 

AOG. 

SUMMARY FOR DIOIDECK BASIN 
. EVENT for; NOVEMBER 
NO. OF storms; 6.78 

RAIN RATE CUMULATIVE 

CUMULATIVE 

4 

+ 

4 

4 

4 

4 

4 

(SECS) 

<M3/S) 


(M/S) 

DISCHARGE(M3> 

RAINFALL(M3) 

4 

4 

4 

4 

900. 

0.00000000 


0.0000000550 

0.0000 

195.6070 

4 

4 

4 

IBOO. 

0.00000000 


0.0000000550 

0.0000 

391.2156 

4 

4 

2700. 

0.00000000 


0.0000001430 

0.0000 

899.7960 

4 

4 

3600. 

0.00000000 


0.0000001430 

0.0000 

1408.3762 

4 

4 

4500. 

0.00000000 


0.0000002310 

0.0000 

2229 . 9290 

4 

4 

5400. 

0.00000000 


0.0000002310 

0.0000 

3051.4819 

4 

4 

6300. 

0.00000000 


0.0000003300 

0.0000 

4225.1289 

4 

4 

7200. 

0.00000000 


0.0000003300 

0.0000 

5398.7759 

4 

4 

BIOO. 

0.00000000 


0.0000004290 

0,0000 

6924.5171 

4 

4 

9000. 

0.00000000 


0.0000004290 

0.0000 

8450.2578 

4 

4 

9900. 

0.00000000 


0.0000005170 

0.0000 

10288.9717 

4 

4 

lOBOO. 

0.00000000 


0.0000005170 

0.0000 

12127.6846 

4 

4 

11700. 

0.00000000 


0.0000006160 

0.0000 

14318.4922 

4 

4 

12600. 

0.00000000 


0.0000006160 

0.0000 

16509.3008 

4 

4 

13500. 

0.00000000 


0.0000008530 

0,0000 

19543.0000 

4 

4 

14400. 

0.00000000 


0.0000000530 

0.0000 

22576.6992 

4 

4 

15300. 

0.00000000 


0.0000005940 

0.0000 

24689.2637 

4 

4 

16200. 

0.00000000 


0,0000005940 

0.0000 

26801.8281 

4 

4 

17100. 

0.00000000 


0.0000005170 

0.0000 

28640.5410 

4 

4 

18000. 

0.00000000 


0.0000005170 

0,0000 

30479.2559 

4 

4 

18900. 

0.00000000 


0.0000004510 

0.0000 

32083.2402 

4 

4 

19800. 

0.00000000 


0.0000004510 

0.0000 

33687.2227 

4 

4 

20700. 

0.00000009 


0.0000003850 

0.0001 

35056.4766 

4 

4 

21600. 

0.00004824 


0.0000003850 

0,0435 

36425.7305 

4 

4 

22500. 

0,00029889 


0.0000003190 

0.3125 

37560.2539 

4 

4 

23400. 

0.00110513 


0.0000003190 

1 . 3071 

38694.7812 

4 

4 

24300. 

0.00229098 


0,0000002420 

3.3690 

39555.4570 

4 

4 

25200. 

0.00381064 


0.0000002420 

6.7986 

40416.1289 

4 

4 

26100. 

0.00495621 


0.0000001760 

11.2592 

41042.0742 

4 

4 

27000. 

0.00576055 


0.0000001760 

16.4437 

41668.0195 

4 

4 

27900. 

0.00598982 


0.0000000990 

21.8345 

42020.1133 

4 

4 

28800. 

0.00598133 


0.0000000990 

27.2177 

42372.2070 

4 

4 

29700. 

0.00595877 


0.0000000330 

32.5806 

42489.5703 

4 

4 

30600. 

0.00591150 


0.0000000330 

37,9009 

42606.9375 

4 


4444444444444444444444444444444444444444444444444444444444444444444+444+4444444444444444444444444444 

++++++++++++++++++++++++++++++++++++++++++++++4++++++++++++++++++++++++++++++++++++++4++++++++++++++ 
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4+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


POLLUTANT LOADING: DIVIDECK BASIN 
EVENTS of: NOVEMBER 
RPC CONCENTRATION DATA 


LAND use: rural- 0.92» LAC= 0.01» LAR= 0.06> MAC- 0.00» HAR= 0.001 CBD= 0.00 

BOD COD 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


45.245 


1.812 


10.522 


TOTAL KJELDAHL 
NITROGEN 


0.210 


TOTAL 

PH0SPDR0U5 


0.107 


COG CONCENTRATION DATA 


pollutant: 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


16.387 


COD 


31.45 


TOTAL KJELDAHL 
NITROGEN 


0.446 


TOTAL TOTAL 

NITROGEN PHOSPOROUS 


0.471 


0.059 


4444444444444444444444444444444444444444444444444444444444444444444+444444444+4444444444444444444444 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
4 


4 
4 
4 
4 
4 

LAND use: SHOP CTR- O.OOf STRP COM- 0.01 J HR RES= 0.001 TH/GA- 0.00> MDSFR- 0.00 4 

LLSFR= 0.06) CON= 0.00) FOR= 0.82! CT AG- 0.00) MT AG- 0.00) PAST- 0.10 4 

4 
4 
4 
4 
4 
4 
4 
4 
4 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


4 

4 

4 

4 

4 

4 

TIME 

niSCHARGE 

RUNOFF 

AUG 

AUG. 

SUMMARY FOR DIUIDECK BASIN 
. EUENT for: DECEMBER 
NO. OF storms: 7.46 

RAIN RATE CUMULATIUE 

CUMULATIUE 

+ 

+ 

+ 

+ 

4 

4 

4 

4 

<SECS) 

<M3/S) 


<M/S) 

DISCHARGE(M3) 

RAINFALL(M3> 

4 

4 

4 

4 

900. 

0.00000000 


0.0000000380 

0.0000 

135.1472 

4 

4 

4 

IBOO. 

0.00000000 


0.0000000380 

0.0000 

270.2944 

4 

4 

2700. 

0.00000000 


0.0000001140 

0.0000 

675.7361 

4 

4 

3A00. 

0.00000000 


0.0000001140 

0.0000 

1081.1777 

4 

4 

4500. 

0.00000000 


0.0000001820 

0.0000 

1728.4618 

4 

4 

5400. 

0.00000000 


0.0000001820 

0.0000 

2375.7458 

4 

4 

6300. 

0.00000000 


0.0000002520 

0.0000 

3271.9854 

4 

4 

7200. 

0.00000000 


0.0000002520 

0.0000 

4168.2246 

4 

4 

8100. 

0.00000000 


0.0000003250 

0.0000 

5324.0889 

4 

4 

9000. 

0.00000000 


0.0000003250 

0.0000 

6479.9536 

4 

4 

9900. 

0.00000000 


0.0000004010 

0.0000 

7906.1123 

4 

4 

10800. 

0.00000000 


0.0000004010 

0.0000 

9332.2715 

4 

4 

11700. 

0.00000000 


0.0000004770 

0.0000 

11028.7246 

4 

4 

12600. 

0.00000000 


0.0000004770 

0.0000 

12725.1777 

4 

4 

13500. 

0.00000000 


0.0000010450 

0.0000 

16441.7266 

4 

4 

14400. 

0.00000000 


0.0000010450 

0.0000 

20158.2754 

4 

4 

15300. 

0.00000000 


0.0000005450 

0.0000 

22096.5703 

4 

4 

16200. 

0.00000000 


0.0000005450 

0.0000 

24034.8672 

4 

4 

17100. 

0.00000000 


0.0000005070 

0.0000 

25838.0156 

4 

4 

18000. 

0.00000000 


0.0000005070 

0.0000 

27641.1641 

4 

4 

18900. 

0.00000000 


0.0000004610 

0.0000 

29280.7129 

4 

4 

19800. 

0.00000000 


0.0000004610 

0.0000 

30920.2617 

4 

4 

20700. 

0.00000000 


0.0000004240 

0.0000 

32428,2207 

4 

4 

21600. 

0.00002987 


0.0000004240 

0.0269 

33936,1797 

4 

4 

22500. 

0.00032761 


0.0000003860 

0,3217 

35308.9922 

4 

4 

23400. 

0.00174515 


0.0000003860 

1.8924 

,36681.8008 

4 

4 

24300. 

0.00527335 


0.0000003400 

6.6384 

37891.0117 

4 

4 

25200. 

0.01220103 


0.0000003400 

17,6193 

39100,2266 

4 

4 

26100. 

0.02257541 


0.0000003030 

37.9372 

40177.8477 

4 

4 

27000. 

0.03679891 


0.0000003030 

71.0562 

41255.4687 

4 

4 

27900. 

0.05299573 


0.0000002650 

118.7523 

42197.9414 

4 

4 

28800. 

0.06991061 


0.0000002650 

181.6719 

43140.4180 

4 

4 

29700. 

0.08294231 


0.0000002190 

256.3200 

43919,2930 

4 

4 

30600. 

0.09053081 


0.0000002190 

337.7977 

44698.1680 

4 

4 

31500. 

0.09084085 


0.0000001820 

419.5545 

45345.4531 

4 

4 

32400. 

0.08571155 


0.0000001820 

496.6949 

45992.7344 

4 

4 

33300. 

0.07610982 


0.0000001440 

565.1938 

46504,8711 

4 

4 

34200. 

0.06476632 


0.0000001440 

623,4835 

47017.0078 

4 

4 

35100. 

0.05345680 


0.0000000980 

671 ,5946 

47365.5469 

4 

4 

36000. 

0.04372070 


0.0000000980 

710.9432 

47714,0820 

4 

4 

36900. 

0.03581090 


0.0000000610 

743.1730 

47931,0273 

4 

4 

37800. 

0.02945849 


0.0000000610 

769,6857 

48147.9766 

4 

4 

38700. 

0.02438080 


0.0000000230 

791 ,6284 

48229.7773 

4 

4 

39600. 

0.02032214 


0.0000000230 

809.9183 

48311.5742 

4 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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4444444444444444444444444444444444444444444444444444444444444t44444444444444444444444444444444444444 
4 4 
4 4 
4 POLLUTANT LOADING! DIVIPECK BASIN 4 
4 EUENTS of; DECEHBER 4 
4 RPC CONCENTRATION DATA 4 
4 4 
4 4 
4 LAND use; RURAL= 0.925 LAC= 0.015 LAR» 0.065 HAC» 0.005 HAR= 0.005 CBD- 0.00 4 


pollutant; 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


TOTAL KJELDAHL 
NITROGEN 


TOTAL 

PH0SP0RQU5 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 


COG CONCENTRATION DATA 


LAND use; shop CTR= 0.005 STRP COM= 0.015 HR RES- 0.005 TH/GA- 0.005 MPSFR* 0.00 

LLSFR- 0.065 CON- 0.005 FOR= 0.825 CT AG- 0.005 MT AG- 0.005 PAST- 0.10 


pollutant; 

LOAD-KGS. 


SUSPENDED 

SOLIDS 


TOTAL KJELDAHL 
NITROGEN 


TOTAL 

NITROGEN 


TOTAL 

PHOSPOROUS 


4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
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+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++4 


44+44+44++++++++++++++++4+44+44++++++4444+44444444+444444444444+44444444444444t4444444444444+4444444 


pollutant: 

LOAO-KGS. 


BOD 


COD 


ANNUAL LOADING - RPC DATA + 

4 
4 

TOTAL 4 

PHOSPOROUS 4 

4 
4 

30,667 4 

4 
4 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

4 
4 
4 


SUSPENDED 

SOLIDS 


129B1 .008 


520.001 3018.782 


TOTAL KJELDAHL 
NITROGEN 


60.272 


ANNUAL LOADING - COG DATA 
COD 


pollutant: suspended cod total KJELDAHL TOTAL TOTAL 4 

LOAD-KGS. SOLIDS NITROGEN NITROGEN PHOSPOROUS 4 

4 
4 

4701.627 9023.32 128.064 135.248 16.879 4 

4 
4 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444444 
44444444444444444444444444444444444444444444444444444444444444444 44444444444444444444444444444444444 



0001 


0002 

0003 

0004 


0005 

0006 
0007 
OOOFJ 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
OOIB 

0019 

0020 
0021 

0023 

0024 

0025 

0026 
0028 

0029 

0030 

0031 

0032 


TABLE AIII: LISTING OF ECHOS, VERSION IP 


DIMENSION RAIN (25) tDURN( 25) r0UTQ(3rl00) y00L(2»100) fUPS(4rl00) » 
1QC(15) fO(lOO) rP(lOO) rHC(lO) fX(lOO) fCTRdS) r IXLN ( 100 ) » OUTQC ( 3 >100 ) 
2>XL0AD(12r 12) >FRAC(12) >T0TL0D(10) >A0G(12) >ARC0G(12) >ARRPC(6) 

3 > CYRLOD < 10 ) r RYRLOD (10) 

DATA ISTARf IPLUS/'*' > '+ V 

DOUDLE PRECISION BASIN r STREAM > BASNAM > DATE 

WRITE (7>411) 

READ HYDROLOGIC PARAMETERS FOR OOERLAND FLOW FROM MASTER FILE 

CALL ASSIGN ( 1 > ' DXl : LINK ♦ DAT ' > 12 ) 

CALL ASSIGN ( 2 r ' DXl : DATA ♦ DAT ' > 12 ) 

CALL ASSIGN ( 3 r ' DXl I RAIN * DAT ' > 12 ) 

CALL ASSIGN ( 5 > ' DXl : CON . DAT ' > 1 1 ) 

PX=0* 

DO 15 IJ=lrlO 
CYRL0D(IJ)=0* 

15 RYRLOD (IJ)=0» 

20 PX^PX+1* 

READ (2>412) NS > D1 > J > T1 > HI 
DO 199 I=1>NS 
READ (1>417) BASIN 

READ (2>413) XL > SL > EN > W > DETEN > CAP > PERC > FI 
READ (3>415) PCT > STORMS >DATE > J1 > TOT 
READ (3 >414) (RAIN(Il) >DURN(I1) >I1=1>J1) 

JXX==J1 + 1 

IF (I .EQ* NS) GO TO 30 
DO 25 MN=1>JXX 
25 BACKSPACE 3 
30 API==CAP)KPCT 

IF (PX *NE* 1*0) GO TO 100 
DO 50 I10=l>100 
50 IXLN(I10)=ISTAR 
WRITE (7>418) IXLN 

WRITE (7>436) BASIN > XL > SL > EN > W f DETEN > CAP > API > PERC > FI 
WRITE (7>435) IXLN 

STATEMENTS 100 TO 199 COMPUTE OUERLAND FLOW 
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0112 


S=S 

-SriELT 



0113 


IF 

(S *LT* 

0. ) 

s=.-o* 

0115 


IF 

(S *GT» 

CAP) 

S=CAP 

0117 


IF 

(K *LT. 

J) GO TO 199 

0119 


IF 

<0 *EQ* 

0* ) 

GO TO 195 

0121 

195 

0UTQ<ItK1)= 



0122 


K^O 




0123 


Kl = 

Kl + 1 



0124 

199 

CONTINUE 




c 

C LINES 200 TO 299 PERFORM CHANNEL ROUTING 
C 


0125 


READ <2f438) BASNAM»AREA 

0126 


TIME=0« 

0127 


QCUM=0* 

0128 


J1C==<T0T/(J)KD1)) 

0129 


READ (2>419) NODES 

0130 


DO 299 M=lr NODES 

0131 


READ (2t 420) XLC r SLC r ENC > WC 

0132 


READ (1»421) KUPSrKOVL, STREAM 

0133 


IF <PX .NE* 1*) GO TO 288 

0135 


WRITE <7r418) IXLN 

0136 


WRITE (7f422) STREAM » XLC f SLC » ENC t WC 

0137 


WRITE (7^435) IXLN 

0138 

288 

DO 293 IJK=lrlOO 

0139 

293 

IXLN<IJK)«IPLUS 

0140 


WRITE (7^418) IXLN 

0141 


WRITE (7M35) IXLN 

0142 


WRITE <7r443) 

0143 


WRITE <7r440) BASNAM » DATE r STORMS 

0144 


WRITE <7»443) 

0145 


WRITE (7>441) 

0146 


WRITE (7^443) 

0147 


IF (KUPS *Ea. 0) GO TO 203 

0149 


DO 202 J2=1>KUPS 

0150 


READ (lr419) K7 

0151 


DO 202 13=1 fJlC 

0152 


UPS(J2fI3) = 0UTQC(K7fI3) 

0153 

202 

CONTINUE 


9e“-V 





0154 

203 

IF (KOOL *EG. 0) GO TO 205 

0156 


DO 204 J2=1>K00L 

0157 


READ (lr419) K7 

0158 


DO 204 14=1 fJlC 

0159 


00L(J2»I4) = 0UTQ(K7rI4) 

0160 

204 

CONTINUE 

0161 

205 

ALFAC = (SLC***5/ENC) 

0162 


TEMPO = Dl>lcJ 

0163 


T=0* 

0164 


DO 210 15 = IflO 

0165 


HC(I5) = 0* 

0166 

210 

QC(I5) = 0* 

FORTRAN 

10 

001C-03A 

0167 


K=0 

0168 


SEGLEN = XLC/10* 

0169 


L3=0 

0170 


K2 = 1 

0171 


DO 299 16 = IrN 

0172 


K=K+1 

0173 


T=T+D1 

0174 


IF (T *GT* TEMPO GO TO 219 

0176 


GO TO 225 

0177 

219 

K2=K2+1 

0178 


TEMPO = TEMPC+(D1*J) 

0179 

225 

DO 228 L=lflO 

0180 


T0TFL0=0* 

0181 


IF (L .NE. 1) GO TO 213 

0183 


IF (KUPS *EQ4 0) GO TO 252 

0185 


DO 206 J6=1tKUPS 

0186 

206 

T0TFi:0=T0TFL0+ ( UPS ( J6 » K2 ) /WC ) 

0187 


GO TO 252 

0188 

213 

T0TFL0=QC(L-1) 

0189 

252 

IF (KOOL .EQ* 0) GO TO 229 

0191 

242 

DO 207 I7=1>K00L 
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0192 207 TOTFLO:=T0TFLO+(OOL(I7fK2) )/(10*>KWC) 

0193 229 IF (HC(L) *LE. 0») GO TO 251 

0195 IF ( (TOTFLO-QC<L) ) *EQ* 0*) GO TO 231 

0197 HC(L) = HC<L) + (D1)K< <TOTFLO““QC(L) )/SEGLEN) ) 

0198 231 IF (HC(L) .LE. 0*) GO TO 281 

0200 241 QC(L) = ( ALFAOI? < HC ( L) **1 . 67 ) ) / ( < 2 . )|cHC ( L) + 1 ♦ ) >K)K ♦ 67 ) 

0201 GO TO 228 

0202 251 IF (TOTFLO *LE» 0*) GO TO 281 

0204 HC<L) = D1)K<T0TFL0/SEGLEN) 

0205 253 IF <HC(L) ♦LE* 0») GO TO 281 

0207 GO TO 241 

0208 281 HC(L)=0» 

0209 0C(L)=0* 

0210 228 CONTINUE 

0211 291 IF (K .LT* J) GO TO 299 

0213 L3=L3+1 

0214 0UTQC<M>L3)=(0C(10)*WC) 

0215 IF <M .NE. NODES) GO TO 299 

0217 DURT0T=0* 

0218 IC0UNT=0 

0219 RNCUM=0* 

0220 TIME=TIME+(D1*J) 

0221 IJK=0 

0222 294 IJK=IJK+1 

0223 IF <IJK *GT* Jl) GO TO 295 

0225 DURTOT=DURTOT+DURN ( I JK ) 

0226 IF (TIME *LE. DURTOT) GO TO 295 

0228 IC0UNT=IC0UNT+1 

0229 GO TO 294 

0230 295 riURTOT=DURTOT“-DURN(IJK) 

0231 IF (ICOUNT *EQ* 0) GO TO 297 

0233 DO 296 IJL=lf ICOUNT 

0234 296 RNCUM==RNCUM+ < AREA*DURN ( I JL) *RAIN ( IJL) ) 
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FORTRAN 10 


001C-03A 


PAGE 


0235 

0236 

0238 

0239 
02-10 

0241 

0242 

0243 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 


297 DIFF=TIME~DURTOT 

IF (DIFF *LT* 0.) DIFF=0, 

ICOUNT-ICOUNT+1 

RNCUM=RNCUM+(AREA*DIFF)KRAIN(ICOUNT) ) 

OCUM=QCUM+ ( D 1 * JJKOUTQC ( M » L3 ) ) 

WRITE (7f442) TIME?0lJTGC<MfL3) fRAIN( ICOUNT) jQCUMrRNCUM 

511 K==0 

IF (16 *NE* N) GO TO 299 
WRITE (7»443^ 

WRITE (7f435) IXLN 
WRITE <7»435) IXLN 
READ <5 t 434) ( ARRPCC ID r IL=1 r 6 ) 

JK=0 

512 JK=JK+1 

IF (JK *EQ* 7) GO TO 515 

READ <5^444) <XLOAD( JK> JL) > JL=1 r5) 

GO TO 512 

515 DO 518 IL6=1>6 
518 FRAC(IL6)=ARRPC(IL6)/AREA 
DO 520 JP=1>5 
TEMP=0» 

DO 520 JUL=1»6 

TEMP=TEMP+FRAC ( JUL ) * XLO AD ( JUL » JP ) 

520 AOG(JP)«TEMP 
DQ 530 IJ=1»5 


0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 


T0TL0D(IJ)=0*001)KST0RMS)KQCUM)KA0G<IJ> 
530 RYRLOD ( IJ ) =RYRLOD ( I J) +T0TL0D < IJ ) 
WRITE (7»418) IXLN 
IXLN 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


<7f435) 

(7r443) 

(7r445) 

(7r443) 

(7^447) 

(7t443) 

<7»448) 

(7t443) 

(7»449) 


BASNAMtDATE 


(FRACdJSl f IJ3=lf6) 


<T0TL0D(IJ4) »IJ4=1»5) 


0275 


WRITE (7>443) 

0276 


WRITE (7r'435) IXLN 

0277 


WRITE (7»435) IXLN 

0278 


WRITE <7r443) 

0279 


WRITE <7»450) 

0280 


WRITE (7^443) 

0281 


READ (5f434) ( ARCOG < ID r IL=1 » 1 1 ) 

0282 


DO 540 IL1=1»11 

0283 

540 

FRAC(IL1)=ARC0G<IL1)/AREA 

0284 


JK=0 

0285 

545 

JK=JK+1 

0286 


IF (JK ♦EQ* 12) 60 TO 550 

0288 


READ (5 >444) ( XLOAD ( JK > JL ) > JL=1 > 5 ) 

0289 


GO TO 545 

0290 

550 

DO 560 JP=1>5 

0291 


TEMP=0* 

0292 


DO 560 JUL^lrll 

FORTRAN 

10 

001C-03A 

0293 


TEMP=TEMP+FRAC( JUL) t XLOAEKJULrJP 

0294 

560 

AOG( JP)=TEMP 

0295 


DO 570 IJ==1>5 

0296 


TOTLOD(IJ) ~0.001*ST0RMS)KaCUM*A0G(: 

0297 

570 

CYRLOD ( IJ ) =CYRLOD ( IJ ) +TOTLOD ( IJ ) 

0298 


WRITE (7>452) <FRAC< IJl ) > IJ1=1 > 11 ) 

0299 


WRITE (7 >443) 

0300 


WRITE (7>453) 

0301 


WRITE (7>443) 

0302 


WRITE (7>454) (T0TL0D(IJ2) >IJ2=1>5 

0303 


WRITE (7>443) 

0304 


WRITE <7>435) IXLN 

0305 


WRITE (7>435) IXLN 

0306 

299 

CONTINUE 

0307 


IF <PX *EQ* 12.) GO TO 300 

0309 


REWIND 1 

0310 


REWIND 2 
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0311 

REUINE 

5 


0312 

GO TO 

20 


0313 

300 WRITE 

(7f418) 

IXLN 

0314 

WRITE 

(7f435) 

IXLN 

0315 

WRITE 

(7f443) 


0316 

WRITE 

(7f439) 


0317 

WRITE 

(7f443) 


0318 

WRITE 

(7f448) 


0319 

WRITE 

(7f443) 


0320 

WRITE 

(7f449) 

(RYRL0D(IJ4) fIJ4=1f5) 

0321 

WRITE 

(7f443) 


0322 

WRITE 

(7f435) 

IXLN 

0323 

WRITE 

(7f435) 

IXLN 

0324 

WRITE 

<7f443) 


0325 

WRITE 

(7f416) 


0326 

WRITE 

(7f443) 


0327 

WRITE 

(7f453) 


0328 

WRITE 

(7f443) 


0329 

WRITE 

(7f454) 

(CYRL0D(IJ2) fIJ2==1f5) 

0330 

WRITE 

(7f443) 


0331 

WRITE 

(7f435) 

IXLN 

0332 

WRITE 

(7f435) 

IXLN 


0333 411 FORMAT (/////» 41X r ' ECOSYSTEMS HYDROLOGIC SIMULATOR (ECHOS) ' r/»57X> 

I'UERSION IP'f//) 

0334 412 FORMAT < I2fF4*0r I3 tF8.0>F7*5) 

0335 413 FORMAT (F8 ♦ 2 ^ 2F5 ♦ 3 » F7 . 0 , F9 ♦ 8 , ,3F1 1 *7) 

0336 414 FORMAT (F10*9»F6.0) 

0337 415 FORMAT (2F5 * 2 » AS r 12 » F7 ♦ 0 ) 

0338 416 FORMAT ( lOX 36X r ' ANNUAL LOADING ~ COG DATA ' , 37X r ' + ' ) 

0339 417 FORMAT (A8) 

0340 418 FORMAT <//» lOXf lOOAl ) 

0341 419 FORMAT (12) 

0342 420 FORMAT ( F8 *2 » 2F5 . 3 » F4 . 2 ) 

0343 421 FORMAT (2I2rA8) 

0344 422 FORMAT ( lOX f IX f ' STREAM SEGMENT: ' f A8f73Xf f/f IOXf f 

1Xf'LENGTH= ' fF8*2f 'METERS' fIOXf 'SL0PE= ' fF5*3f IOXf 'N= 'f 
2F5*3f10Xf 'WIDTH= ' fF4*2f 'METERS' f8Xf 't'> 
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FORTRAN IV 


V01C-03A 


PAGE 007 


0345 

034^> 

0347 


034B 

0349 

0350 


0351 


0352 


0353 

0354 

0355 


0356 

0357 

0358 


0359 

0360 

0361 

0362 


434 FORMAT <6F10.1) 

435 FORMAT (IOXtIOOAD 

436 FORMAT ( lOX IX r " SUBBASIN ; ' > AS r 79X r ' >K ' » / r lOX y ' )K ' y IX f ' LENGTH^ ' 

lyF8.2y' METERS' f2Xy ^SLOPE= ' r F5 . 3 y 2X y ' N= ' y F5 ♦ 3 y 2X y ' UIDTH= ' y F7 ♦ 0 y 
2 ' METERS ' y2Xy 'DETENTION^' yF7*5y ' METERS ' y 2X y ' ' y / y lOX y ' ' y IX y 'SOIL 
3 MOISTURE CAPACITY = 'yFll»9y' METERS ' y 5X y ' ANTECEDENT SOIL MOISTUR 
4E== 'yFll*9y' METERS' y4Xy '*' y/ylOXy ')K' rlXy 'PERCOLATION RATE= 'yFll* 
59y' M/SEC' yl4Xy 'FINAL INFILTRATION RATE= 'yFll*9y' M/SEC ' y 6X y '*' ) 

438 FORMAT (A8yF10*l) 

439 FORMAT ( lOXy '+' y36Xy 'ANNUAL LOADING - RPC DATA ' y 37X y ' + ' ) 

440 FORMAT (lOX y '+' y 33X y ' RUNOFF SUMMARY FOR 'yA8y' BASIN'y 
132Xy '+' y/y lOXy '+' y37Xy 'AUG* EUENT FOR: ' y AB y 37X y ' + ' y / y 
210Xy '+' y36Xy 'AVG* NO* OF STORMS: ' y F5 ♦ 2 y 37X y ' + ' ) 

441 FORMAT < lOX y '+' y 9X y ' TIME ' y lOX y ' DISCHARGE ' y 9X y ' RAIN RATE'ylO 
IX y 'CUMULATIUE' y9Xy ' CUMULATIUE ' y 9X y '+' y/y lOXy '+' yBXy ' <SECS) 

2' y llXy ' (M3/S) ' y 12Xy ' (M/S) ' y 12Xy ' DISCHARGE ( M3 ) ' y6Xy 'RAINFALL 
3(M3) ' y7Xy '+' ) 

442 FORMAT (lOXy ' + ' y 7X y F7 ♦ 0 y 6X y F12 ♦ 8 y 8X y F12 ♦ 10 y 4X y F13 * 4 y 7X y F12 * 4 y 
llOXy '+' > 

443 FORMAT ( lOX y ' + ' y 98X y ' + ' y / y lOX y ' + ' y 98X y ' + ' ) 

444 FORMAT (5F7*2) 

445 FORMAT ( lOX y '+' y 33X y ' POLLUTANT LOADING: 'yA8y' BASIN'y 
132Xy '+' y/y lOXy '+' y40Xy 'EVENTS OF: ' y A8 y 39X y ' + ' y / y lOX y 
2'+' y38Xy 'RPC CONCENTRATION DATA ' y 38X y ' + ' ) 

446 FORMAT (A4) 

447 FORMAT ( lOX y ' + ' y 1 IX y ' LAND USE: RURAL= 'yF5*2y'y LAC= 'yF5.2y'J 

1 LAR= 'yF5*2y'y HAC= 'yF5*2y'y HAR== ' y F5 * 2 y ' y CBD= ' yF5*2y5Xy ' + ' ) 

448 FORMAT (lOXy '+'y9Xy 'POLLUTANT: ' y 7X y ' SUSPENDED ' y 7X y 'BOD' y 
27Xy 'COD' y7Xy 'TOTAL K JELDAHL ' y 7X y 'TOTAL' ylOXy ' + ' y/y lOXy 
3'+' y9Xy 'LOAD-KGS* ' y8Xy ' SOLIDS ' y 30X y ' NITROGEN ' y 13X y ' 

4PH0SP0R0US ' y 5X y ' + ' ) 

449 FORMAT ( lOXy '+' y24XyF10*3y4XyF8*2y3XyF8*2y5XyF12*3y9Xy 
2F12*3y3Xy '+' ) 

450 FORMAT ( lOX y ' + ' y 38X y ' COG CONCENTRATION DATA ' y 38X y ' + ' ) 

451 FORMAT (F10*2) 

452 FORMAT < lOX y ' + ' y 9X y ' LAND USE: SHOP CTR= 'yF5*2y'y STRP COM= 'y 
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0363 


0364 

0365 

0366 
''C 


HR RES= 'fF5*2f'f TH/GA= 'fF5»2»'f z.’ SFR= '»F5«2» 
24X>'f'>/>10Xr'+'fl9X»'LLSFR= 'fF5*2>'? CON '»F5*2f'? FOR= 
3F5*2f'r CT AG= 'fF5*2r'i MT AG= 'tF5»2>'? ' rF5 *2f2Xf ' i' ) 

453 FORMAT <10X» ' + ' >9Xf 'POLLUTANT: ' » 7X r ' SUSPENDL"! ' , 7X » 'COD' r7X» 

1' TOTAL KJELDAHL' r7Xf 'TOTAL' ^7Xr 'TOTAL ' r 8X f ' -f '»/» lOX f ' + ' r 
29 Xt 'LOAD-KGS* ' reXr 'SOLIDS' t 20X» 'NITROGEN ' r 13X f ' NITROGEN '» 4Xr 
3'PH0SP0R0US'F3Xr '+') 

454 FORMAT <10X> ' + ' y 24X r FIO * 3 » 4X » F8 ♦ 2 r 2 ( 9X » F8 ♦ 2 ) . ‘ X y FIO . 3 y 3X y ' + ' ) 
999 STOP 

END 
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Milton H. Miller, Chairman 
C Bowie Rose, Sr, Vice Chairman 
Frederick L Dewberry, 

Executive Director 

701 St. Paul Street 
Baltimore, Maryland 21202 
(301) 383-5838 



April 20, 1978 


Mr, Harry Loats 

Technical Director 

ECO Systems International, Inc. 

P. 0. Box 225 

Gambrills, Maryland 21054 

Re: Performance and Application 

of the Magothy (Anne Arundel 
County, Maryland) Hydrologic 
Planning Analysis 

Dear Harry: 

The Baltimore Regional Planning Council's (BRPC) Section 
208 Water Quality Management Planning Program would like to 
sincerely thank ECO Systems International for its technical 
support in hydrologic planning of the Magothy River Basin 
(Anne Arundel County, Maryland). Chosen by BRPC as one of 
three areas of intensive study, the Magothy Basin represents 
one of numerous tidal watersheds along the Chesapeake Bay's 
Western shore. ECO Systems has provided land cover classifi- 
cation and hydrologic planning tools - 

(1) defining spatial location and extent of seven 
land cover types by watershed, 

(2) simulating storm-related runoff, and 

(3) quantifying non-point or diffuse sources of 
water entrained pollution flux, - 

through practical application. 

As the primary goal of Public Law 92-500 is to provide 
"fishable and swimmable waters" in the nation’s waters by 1983, 
the BRPC is presently engaged in the water quality planning pro- 
cess for seven major river basins. Briefly, the process in- 
cludes: (1) identification and quantification of pollutant 

sources for both existing and projected future conditions, 

(2) design of problem abatement alternatives, (3) testing and 
evaluation of alternatives for selection of the most cost- 
effective and feasible one by elected decision-makers. 


Baltimore City Anne Arundel County Baltimore County Carroll County Harford County Howard County State of Maryland 



B-2 


LETTER - Mr, Harry Loats 


- 2 - 


April 20, 1978 


The ECO Systems analysis and resulting methodology allows BRPC to 
evaluate non-point source pollutant flux temporarily and spatially under both 
existing and projected land cover conditions. Sub-basins will be evaluated 
and ranked according to established selection criteria from highest to lowest 
flux generation , This focuses the reduction of pollutants in sub-basins with 
the most critical problems. 

To follow the initial problem analysis for existing conditions, the 
water quality planning process will be composed of the following steps: 

(1) Project land cover changes and compute the resulting pollutant 
load generation and delivery for a selected number of alterna- 
tives 

(2) Establish the impact on the estuarine receiving water quality 
from land generated pollutant flux by application of steady 
state modeling 

(a) existing conditions 

(b) projected future conditions 

(3) Apply broad NPS pollutant control strategies by existing and 
projected pollutant source loads 

(4) Select the most cost-effective control strategy and critical 
areas of primary application. 

In the proposed continuing planning process, ECO Systems hydrologic 
planning methodologies - as applied in the Magothy River - should be used 
throughout the BRPC region for non-point source pollutant quantification. 

The present planning program allows only portions of three river basins to 
be studied with sufficient resources to propose detailed control strategies, 

ECO Systems methodologies require low resource expenditures with most hydro- 
logic data available from existing documentation. Therefore, water quality 
planning staff should be able to apply the methodologies with minimal training. 

Sincerely 


Samuel R. Martin 
Non-Point Source Engineer/ 
Planner 

Baltimore Water Quality 
Planning Program 
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May 1, 1978 


Mr. Harry Loats 
Technical Director 
ECO Systems International, Inc. 
Box 225 

Gambrills, Maryland 21054 


Milton H. Miller, Chairman 
C. Bowie Rose, Sr,, Vice Chairman 
Frederick L Dewberry, 

Executive Director 

701 St. Paul Street 
Baltimore, Maryland 21202 
(301) 383-5838 


Dear Harry: 

I wish to express my appreciation for the work which you 
did on the Magothy River for the Baltimore Region 208 plan- 
ning effort. We received two major benefits from your effort. 
First, the Magothy work significantly advanced the use of 
satellite data as a principal methodology for use by the 
Regional Planning Council in future planning projects and also, 
the water quality and pollutant loading analysis was of ma- 
terial benefit in helping us to understand the mechanics of 
non-point source and estuarine water quality problems. 


Thank you for the special effort which you undertook to 
help us. 


Sincerely yours 



H. Clayton Ervine 
Technical Coordinator 


Baltimore City Anne Arundel County Baltimore County Carroll County Harford County Howard County State of Maryland 
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